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ON THE WIND GENERATION OF WATER WAVES 


OLIVER BUHLERt, JALAL SHATAH*, SAMUEL WALSH, AND CHONGCHUN ZENG§ 


Abstract. In this work, we consider the mathematical theory of wind generated water waves. 
This entails determining the stability properties of the family of laminar flow solutions to the 
two-phase interface Euler equation. We present a rigorous derivation of the linearized evolution 
equations about an arbitrary steady solution, and, using this, we give a complete proof of the 
instability criterion of Miles [15]. Our analysis is valid even in the presence of surface tension 
and a vortex sheet (discontinuity in the tangential velocity across the air—sea interface). We are 
thus able to give a unified equation connecting the Kelvin-Helmholtz and quasi-laminar models 
of wave generation. 


1. Introduction 


In this paper, we seek to address the extremely classical problem of determining how wind 
blowing over the ocean generates waves. Specifically, our main objective is to give a mathematically 
rigorous answer to the question: What must the velocity profile of the wind be in order to give rise 
to persistent waves in quiescent water? 

For our purposes, the air-sea system is modeled by the two-dimensional incompressible interface 
Euler problem. That is, we consider the ocean and atmosphere as immiscible fluids, each evolving 
according to the free surface Euler equations. The assumption of incompressibility is widely adopted 
and reasonable because the Mach number for typical flows is quite small (see, for example, [14]). 
At time t > 0, the air inhabits the region , and the water occupies 17)”. The ocean is finite depth 
with a rigid flat bed. We also follow the common practice of taking the atmosphere region to lie 
below a rigid flat lid; this is justified on the grounds that, if the flow is evanescent at high altitudes, 
its behavior there does not strongly affect the dynamics near the ocean (see, e.g., [14]). Letting 
St := dilf n i9S7(” denote the air-sea interface, we assume moreover that 

n+U^tUSt =:ntUSt=Tx {-h-,h+). 


Here T := ]R/27rZ is the one-dimensional torus. The constants h± > 0 are fixed and describe the 
location of the ocean bed and atmospheric lid. For simplicity, we will suppose that St is given as 
the graph of a smooth function rj = ri{t,xi). This is not strictly necessary, but simplifies many of 
the computations. 

The velocity field v = v{t, x) € and pressure p = p{t, x) £ M satisfy 


(1.1a) 

(1.1b) 


Dju -f- - Vn -I- 002 = 0 

in 17, 

o 

II 

> 

in 17, 
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where Dt := dt + is the material derivative, p is the density, and g is the gravitational constant. 
Implicit above is the assumption that p is constant in each fluid region. We point out that the 
momentum equation (1.1a) does not include any turbulent effects — this is a modeling choice that 
we discuss below. 

The motion of the interface is dictated by the kinematic condition which expresses the fact that 
St is a material line 

(1.1c) dt + v±-y is tangent to \ x £ St}. 

Likewise, the rigidity of the ocean floor and atmospheric lid means that 
(l.ld) Wi • 02 = 0 on X 2 = ih±. 

Note that we are using the convention that, for a quantity / with domain fit, f± := /|f 2 ±- Lastly, 
we impose the dynamic boundary condition 

(l.le) p+(t,x) — p-(t,x) = aK(t,x) = —for all x G S't. 

Here cr > 0 is a (fixed) material constant, and /t(t, x) is the signed curvature of St at x. 

The local well-posedness theory for (1.1) has been studied intensively. Naturally, the irrotational 
problem has enjoyed the most attention, but as will become clear, vorticity in the atmosphere 
plays a central role in wind-wave generation. Thus we will limit our discussion to the literature 
concerning the rotational case. With or without vorticity, when surface tension is neglected (i.e., 
taking a = 0), the linearized system is ill-posed (cf., [4]). When cr > 0, the full nonlinear problem 
is locally well-posed (cf., [6, 20]) 

The interface Euler equations possess a large class of nontrivial exact solutions: observe that any 
pair [v, St) with the ansatz 

v{t,x) = [U{x 2 ), Q), = T X {0} 

comprises a time-independent solution of (1.1). These are called laminar or shear flows, and serve 
as a model for the state of the air-sea system before water waves have formed. If the background 
flow is stable, then perturbations will remain in a neighborhood of the equilibrium, meaning that 
the interface does not leave its quiescent state. On the other hand, if it is unstable, then the free 
surface will become deformed at finite amplitude, i.e., persistent surface waves will be born. In this 
way, the study of wind generation of water waves is equated with diagnosing the stability/instability 
of the laminar flows subject to Eulerian dynamics. 

A natural starting point is to consider the situation where the ocean is at rest and the velocity 
of the air is uniform, i.e., U- = 0, U+ = Uq, for some constant Uq. This is the classical Kelvin- 
Helmholtz model, and it is indeed (linearly) unstable whenever cr = 0 and Uq 0 (in a sense that 
we will make precise later). With surface tension, the flow is (linearly) stable for Uq satisfying the 
inequality 

U^<2P±±f^[ga{p.-p+)]^/^. 

P+P- 

(See, e.g., [9]). Setting p± and cr to their physical values, this predicts that the onset of instability 
occurs when Uq > 6.6 m/s, which is an order of magnitude larger than observation suggest. Worse 
still, this instability is first felt at very small wavelength, roughly 0.0I7m; to excite a wave with 
a more typical wavelength of less than 100m, say, requires Uq to be another order of magnitude 
larger. 
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One is forced to admit, therefore, that the Kelvin-Helmholtz instability fails miserably as a model 
of the wind generation of ocean waves. In particular, it is missing some destabilizing mechanism 
inherent in the physical system. Kelvin himself observed this fact in his original article on the topic 
in 1871 (cf. [23]). Since then, the task of finding a suitable replacement has been a fundamental 
problem in geophysical fluid dynamics. The next century saw a succession of competing models (cf., 
e.g., [14] for a summary), but now the majority opinion has largely settled on the quasi-laminar 
model. Put forward by Miles in a series of papers (cf. [15, 16, 17]), the quasi-laminar model views 
the wind generation process as a resonance phenomenon. Briefly, the idea is the following. Since 
p+lP- ^ Ij the atmosphere may be viewed as a perturbation of vacuum. For infinite depth gravity 
water water beneath vacuum, one has the dispersion relation c = \/g/k, where c is the wave speed 
and k is the wavenumber. That is, the formally linearized problem has an exact solution for which 
St is given as the graph of a function proportional to -<=*)_ This suggests that the dispersion 
relation with an atmosphere may have the asymptotic form 

(1.2) c = -f CiC-b O(e^), e 

In view of (1.2), the situation where a critical layer exists is of special interest: in this case the 
unperturbed phase speed \fgjk lies in the range of flow speeds U+ and this may enable a certain 
resonance between the shear flow in the atmosphere and the gravity waves in the ocean, which then 
manifests itself as a linear instability. Formally linearizing (1.1) and searching for a growing mode 
solution of the form v = Miles concluded that (linear) instability occurs if a critical 

layer location exists where y/g/k = [/(a;^) and U"{x 2 ) < 0. This computation was done with 
(7 = 0, meaning that the Cauchy problem for the full system is ill-posed. Nonetheless, if one fixes a 
physically realistic value for the wavenumber fc, one can still study the linear evolution. Doing this, 
Miles’s scheme gives a way to estimate the corresponding growth rate. Most importantly, in contrast 
to Kelvin-Helmholtz, the quasi-laminar model allows for wave generation without unreasonably 
large wind speed. 

Though it is now widely accepted, the quasi-laminar model has been criticized in the applied 
literature for consigning turbulence to a relatively minor role (cf., e.g., the discussion in [14, Chapter 
4]). Indeed, turbulence is absent in the basic equations (1.1), and so its influence is felt only 
implicitly. For instance. Miles uses a logarithmic wind profile in his growth rate computations, 
which is predicted by the theory of turbulent laminar flow over a flat plate. Also, the presence 
of a turbulent boundary layers is used to justify his assumption that C7+(0) = t/_(0) = 0. Since 
the actual air flow near the ocean surface is known to be highly turbulent, many authors have 
argued that the effects of turbulence must be included explicitly in the dynamics of the air flow. 
Nonetheless, both in field observations and laboratory studies [18, 11], the quasi-laminar model has 
proven to be a rather good predictor of wind energy transfer in many settings. Of course, when 
considering these studies, one must take into account the extreme difficulty in obtaining accurate 
readings; even the task of deducing the wind profile is highly nontrivial. In short, there is a great 
deal of uncertainty, but the ideas of Miles have been more-or-less borne out by the subsequent 
half-century of experimentation. 

Summary of results. Having established the background, let us now enumerate the mathematical 
contents of the present work. It consists of two parts. The first is a careful derivation of the 
linearized interface Euler system (1.1) about an arbitrary steady solution, including those with 
a vortex sheet and taking into account surface tension. The second is a rigorous analysis of the 
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behavior of solutions of the Rayleigh stability equation when the background flow has a critical 
layer occurring at an non-inflection point. We elaborate on both of these points below. 

Rigorous linearization. All prior derivations of the linearized equations for the two phase Euler flow 
are done heuristically and in a setting that does not allow surface tension and/or a vortex sheet if 
the background velocity is non-uniform. In [14], for instance, Janssen suggests that the two-fluid 
problem be imagined as a limit of single-fluid problems with smooth heterogeneous density. The 
advantage of this approach is that it obviates the need to fix the domain, since the one phase fluid 
simply occupies the entire plane. However, the limiting process as the density becomes discontinuous 
over the interface is not straightforward. It is hard, for instance, to see even at a formal level how 
the equations governing the dynamics of the free boundary arise when one allows for a vortex sheet 
and surface tension. Moreover, rigorously establishing that the limit of the one phase problems 
converge to a solution of the two phase problem is a very difficult task (cf., e.g., [13] where it is 
carried out in a very specific physical regime). 

Another typical approach is to linearize the Euler equation (1.1a) separately in each fluid region, 
then formally linearize the equation of boundary motion (1.1c). This is precisely what Miles does 
in [15], where again we note that he assumes the continuity of the tangential velocity across the 
interface and works without surface tension. A flow is then considered linearly unstable if there 
exists a growing mode solution of the resulting problem. 

There is one feature that this process lacks to be mathematically rigorous. Notice that we are 
considering a free boundary problem, so the meaning of “linearized operator” is somewhat subtle. 
From a dynamical systems point of view, the only suitable definition of the linearized propagator is 
found by taking the Frechet derivative or Gateaux variation. For this to make sense, the linearized 
problem must first be formulated in a specific function space. In the present setting, this issue 
manifests itself when one takes variations of the velocity field. Performing the standard formal 
linearization of the Euler equation, we can see immediately that the linearized velocity will have 
the wrong boundary conditions on the interface and hence is not in the correct function space. This 
is expanded upon in Section 2; specifically (2.7) and Remark 2.4. 

Our linearization based on the Hodge decomposition leads us to a system consisting of an ODE, 
Rayleigh’s equation, coupled to an algebraic equation related to the dispersion relation. On the 
one hand, in the absence of surface tension and a vortex sheet, this coincides with the one derived 
formally by Miles. Thus our results mathematically corroborate and generalize his. On the other 
hand, taking the background flow to be uniform in the air and water regions, we recover the classical 
Kelvin-Helmholtz instability criteria (see Example 3.1). This means that we can consider simulta¬ 
neously the destabilizing effects of a vortex sheet and a critical layer, and we can see quantitatively 
the individual influence of each one on the leading-order terms of the unstable eigenvalue. 

A proof of Miles’s criterion and generalizations. Using our linearized problem, we look for an 
unstable eigenvalue. At this point, in [15], Miles proceeds with the heuristic argument outlined 
above, concluding that a linear instability will manifest at wavenumber k provided that \fgjk is in 
the range of the wind profile U+, and that it occurs at an altitude where U'f_ is negative. 

Our presentation is the first rigorous treatment of the system studied by Miles and its general¬ 
ization that we have derived. Mathematically, the analysis involves a detailed examination of the 
Rayleigh equation in the presence of an imaginary parameter — the complex wave speed c — that 
couples the ODE to the dispersion relation. 

The Rayleigh equation (also called the inviscid Orr-Sommerfeld equation) arises when investi¬ 
gating the stability of laminar flows in a rigid channel and has an extensive literature. However, 
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the challenges we confront here are quite distinct from those typically encountered. For shear flows 
in a channel, the boundary conditions are simply homogeneous Dirichlet, whereas the free surface 
in our system leads to an inhomogeneous condition on the interface (see also [12]). More signifi¬ 
cantly, for channel flows, most rigorous sufficient conditions for instability are based on bifurcation 
analysis with c near an inflection point of U. In such cases, treating the wave number as the bi¬ 
furcation parameter, one finds that the leading-order part of the Rayleigh equation is in fact not 
singular. The situation is almost the complete opposite in our work: the critical layer instability 
occurs precisely when c is close to some U{x 2 ) with {/"(xj) 0. This means that the crucial part 

of our analysis is near the singularity of the Rayleigh equation where the solutions will develop a 
logarithmic singularity in the limit Re c —^ 0. One may find plenty of asymptotic expansions of the 
Rayleigh equation solutions near such a singularity in the literature, some of which can be rigorously 
justified provided that U is analytic [25]. By contrast, our approach based on more modern tools 
from dynamical systems, requires only that U+ € and is transparently rigorous. Along the way, 
we obtain a result of independent interest characterizing the limiting behavior of solutions to the 
Rayleigh equation as the wave speed c converges to the real axis; see Proposition 4.1. 

This analysis is then used to prove the main contribution found in Theorem 4.1. For a fixed 
wave number k, let Ck denote the wave speed for the corresponding (linear) capillary-gravity wave 
in finite-depth water beneath vacuum (cf. (3.7)). We obtain a sufficient condition guaranteeing the 
existence of an unstable wave speed lying in a neighborhood of Ck that is valid even with a vortex 
sheet and incorporating surface tension. We are, moreover, able to determine the leading-order 
terms of the unstable eigenvalue, which gives a means of estimating the energy transfer rate; see 
(4.30) and (4.31). Additionally, in Proposition 3.1, we obtain an exact (though implicit) formula 
for the dispersion relation analogous to (1.2). All of these results agree with the formal analysis of 
Miles when the background velocity is continuous over the interface and surface tension is neglected. 
Moreover, in Lemma 3.1, we show that if a sequence of the unstable wave speed Ck^e converges to 
Cfe as e = — —> 0, then Ck must be in the in the range of 17+. Roughly speaking, this says that a 
critical layer is necessary for the generation of instability. 

It is also worth comparing the linear instability in the channel flow problem with the water-air 
interface problem. In the case of channel flow between rigid boundaries, where the wave number is 
often used as the parameter, the instability occurs at isolated wave numbers and thus the linearly 
unstable waves are superpositions of plane waves of isolated wave numbers. On the contrary, for 
the water-air interface problem considered here, the instability due to critical layers occurs at all 
wave number k in certain intervals (or a union of intervals). In particular, with surface tension, 
sufficiently large or small wave numbers are always stable. Without surface tension, sufficiently 
long waves are always stable, while the instability of short wave is affected by both the critical 
layers and Kelvin-Helmholtz instability. The linearly unstable waves are superpositions of waves 
with wave numbers ranging in an interval. 


Unstable waves with critical layer at an inflection point. The critical layer analysis is valid under 
the assumption that the shear profile U is of class (7^ in the air region. As an application of the 
exact dispersion relation, in Section 5 we construct an explicit piecewise linear U (thus U" is a 
(5-function) that are linearly unstable. However, it is shown that, for a certain range of parameters, 
Ck is in the range of C/, but occurs away from the mass of the ^-function, hence U" vanishes there. 
This implies that the instability is not arising from resonance with the critical layer but instead 
from the immensity of l?7"j. For such waves, we have that wave speed c is at an 0(1) distance from 
Ck, whereas the expansion (1.2) is only valid up to 0(e). In fact, demonstrates the sharpness of our 
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necessary condition for instability in Lemma 3.1: there we require that 17+ € C^, whereas these 
profiles are of class 


2. Linearization 

In this section, we will derive the linearization of the Euler equation (1.1) about a steady solution 
with graph geometry. That is, we consider a solutions of the form 

(2.1) Sf = S = {X2 = v = v*{xi,X2), 

for some smooth wave profile rf. One should think of {v*,S) as representing a traveling wave 
observed from a moving reference frame so that it appears stationary. For the quasi-laminar model, 
we are specifically interested in shear flows, i.e. solutions of the form 

(2.2) 5" = {a;2 = 0}, v*{xi,X2) = {U{x2),0), xi gT, X2 & [—h-,h+] 

where 

U{X2) = U-{X2)X{X2<0} + U+{X2)X{X2>0} 

and U± are smooth functions on ± 0:2 G [0,/i±]. Note that we are allowing v* to have a jump 
discontinuity over S. Eventually we will consider those shear flows with 17_ = 0, i.e. the water will 
be assumed to be stationary while there is wind in the air. For those shear flows, the pressure is 
hydrostatic, 

(2.3) Vpl = -gp±e2. 

It is elementary to confirm that all shear flows are solutions of the Euler system. The existence of 
traveling waves where S is not flat has been established by many authors in various regimes (cf., 
e.g., [1, 24]). 

Admissible spaces and orthogonal decompositions. Before we begin, we must introduce the 
spaces in which the problem is formulated. There are several results that prove the local well- 
posedness of the Cauchy problem for the interface Euler equation (1.1) (cf. [6, 20]). These consider 
smooth velocity fields and surface profiles, for example v G for s > 5/2. To establish linear 

instability, we will assume that the background flow is smooth and seek solutions of the linearized 
problem that grow in the norm. 

With that in mind, for any time t > 0, we consider velocity fields v belonging to the Sj-dependent 
space 

X(St) := {v: L^(T X [-h., h+],R^; p dx) ] V • i;± = 0 in Iff, 

• 02 = 0 on X 2 = ±/i±, v+ ■ N+ + V- ■ N- = 0 on St} 

where N± are the unit outward normals for nf. Note that they obviously satisfy A+ + A_ = 0. 
The boundary conditions on St included in the definition of X(S't) are meant to guarantee that 
V • V = 0 in the distribution sense on T x [—h+j. For a divergence free vector field v in L'^, 
its normal component v ■ N is well-defined in H~^{St), as St is smooth (cf., e.g., [22]). X[St) is a 
subspace of L^(T x [—/i-, h+\,pdx) and its orthogonal complement X(5't)-*- is given by 

^{St)^ = {v = -Vq 1 q = q+Xn+ + 9-Xny > 9± e and p+q = p-q on 

This can be seen from the Hodge decomposition as described in [19, 20]. For any X G L^(T x 
[-h-,h+],pdx), 

(2.4) there exists w G 'K{St) such that X = w — Vq and — Vq G X{St)^ 
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where q is determined (uniquely up to a constant) by 
r-Ag = V-X (Tx (-h_,h+))\5t 

(2.5) I q±\s, = := ■ {X+ - VA;V • A+) + A_ • (A_ - VAlV • A_)) 

[V(7± • e2\x2=±h± = X± ■ 02 

Here we write (A±)“^ to denote the inverse Laplacian on with zero Dirichlet boundary condition 

on St and zero Neumann boundary condition on {x 2 = ±h±}. Likewise, we let 

AT := —A4 + —AT. 

P+ P- 

with A± being the Dirichlet-to-Neumann operator on St associated to A± with zero Neumann data 
on {x 2 = ±h±}. Again, we note that for X as above, 

N± ■ (A± - VA±V • A±) e H-^St), 

since X± — VA^^V • X± is in and divergence free while St is assumed to be smooth. 


Linearized evolution equation. To compute the linearized equations, let (5t(a), u(a, t, x)) along 
with the pressure p(a, t, x), be a one-parameter family of solutions of the Euler equation (1.1) that 
coincide with the steady state (S', v*) at a = 0. Since St is close to S (the graph of r]*) for small a, 
we may represent St as the graph of a unique function rj with r]\a=o = V*- 

St = {X2 = r]{a,t,Xi), xi € R}. 


Eventually we will specialize this to the case of a shear flow that is periodic in xi, but for now we 
work in the general setting. 

Erom the boundary condition (l.le), we have (suppressing the dependence on t and a) 


P+ {xi,v{xi)) 


- p-{xi,r]{xi)) = aK = 


(1 -b 772^3/2- 


Let z{t,Xi) := {dati)\a=o, which is the component of the linearized solution corresponding to the 
variation of the interface. Differentiating the above equality we obtain 


( 2 . 6 ) 


aK'{-q*)z = daP+{t,Xi,p*{xi)) - daP-{t,Xi,T]*{xi)) -f {da;2P*+ “ da:2P-){xi, T]* {xi)) Z 


where 








Since v{a,t) G X(St(a)), it follows that 


(v+{xi,r]{xi)) - v-{xi,r]{xi))'j ■ = 0. 

In the above equation, we are again suppressing the dependence on a and t. 
and evaluating at a = 0, we obtain 

(2.7) {ldavl + ldx 2 V*lz) ■ {-dxiP*,lV = lx* ■ eijdxiZ on S. 


Differentiating in a 


Here, for a function / defined on T x [—we write |/] := (/+ — /-)|s. Observe that this 
computation shows that in general the linearized velocity field daV is not in X(S') even though daV± 
is divergence free in and daV± • 02 = 0 along {x 2 = ±h±}; see also Remark 2.4. Our next step 
is therefore to decompose daV into two components, one lying in X(S') and the other in X(S')-’-: 

do,v = Y + Vr, r G X(S') and vr G X(S')^. 





O. BUHLER, J. SHATAH, S. WALSH, AND C. ZENG 


From (2.4), (2.5), and (2.7) 


( 2 . 8 ) 


AF = 0 in rj 

dxj'± = 0 on {x 2 = ±h±}, 

1 


^v*lz ■ - lvllda;^z) on S, 

where Q is the non-normalized Dirichlet-Neumann operator for S' in T x [—/i_, 


r± = —g-^ 

p± 


e±:=tA+^^A4, g-.= —g+ + —g- 

VI p+ 

Likewise, linearizing the momentum equation (1.1) we find that daV = Y + VF satisfies 
(2.9) Ft+ (?;*• V)r + (i;* • V)VF + (y • V)z;* + (VF • V)?;* + VP = 0 in fl, 
where Y = (^ 1 ,^ 2 )^ and 

P = —Pa + Tj. 

P 

Since Y G X(S), taking the divergence and the normal component along X 2 = ±h± of (2.9) and 
using (2.6) and (2.8), we can determine P by solving 


{ -AP = V • ((i;* • V)y + {v* ■ V)vr + (Y ■ W)v* + (VF • V)?;*) in fl 
dx^P = 0 on {x 2 = ±/i±} 

P+P+ - p-P- = -{dx^pX - dx 2 P-)z + (JK'{r]*)z on S 

Finally, because St is given as the graph of 77 , the fact that the normal velocity of the fluid 
interface coincides with the normal component of the velocity field along the interface translates to 
the following statement 

77t(xi) = v{xi,-q{xi)) ■ (-77'(a;i),l)^. 

Linearizing this equality gives 

(2.11) Zt + {v±-ei)dx^z={do,v± + dx 2 V±z)-{-dx:^'q*,lY' on S. 

Due to (2.7), the above equation does not depend on the choice of + or — sign. 

Evolution equations (2.9) and (2.11) along with (2.8) and (2.10) form the linearized system of 
the two phase fluid Euler equation at an arbitrary steady solution {v *, S) with graph geometry. 

In the present work, we are mainly interested in the stability of shear flows under periodic 
perturbations. It is therefore useful to record how these equations simplify for such flows. If we 
take {v*,S) to be given as in (2.2), then in particular S' = T x {0}, so G = Af and k'( 0) = 3^^. 
Thus (2.9) and (2.II) become 

(2.12) + +P'y 2 ei + PVF,,,+r,,,P'ei+VP = 0 in S 

(2.13) Zt + U±{0)zxi=daV±-e2 = Y 2 + dx 2 ^± on {a ;2 = 0}, 


where F and P are determined from 

, , , AF = 0 in n 

(2.14) 

and 


1x2 —— 0 


F±|s = ^(P-(0)-P+(0))S,,Ar-iz, 


f-AP = 2U'{dx,Y2 + Tx,x2) in ft 

< dx^P = 0 on {x 2 = ±h±} 

[p+P+ - P-P- = 9iP+ - P-)z + crzxtxt 


(2.15) 


on S = {x 2 = 0}. 
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Eigenvalues and eigenfunctions. Notice that the coefficients in the linearized system (2.12), 
(2.13), (2.14), and (2.15) depend only on X 2 - Therefore, each Fourier mode is decoupled from 
the other modes. Consider solutions of the linearized system taking the form 


(2.16) 


(z,r,r,p) = k e z\{o}. 


which represents and eigenfunction for the linearized system corresponding to an eigenvalue —ike. 
Clearly, the existence of a solution of this type with Im c > 0 immediately implies exponential linear 
instability. 

In what follows, we derive a linear system for the above unknowns, fixing the Fourier mode k. 
Before doing that, let us record the symbol for the Dirichlet-Neumann operators that we employ: 

Af±{k) = \k\ tanh {\k\h±) 

(2.17) ^ 1 _ 1 

—Af+{k) + - 
P+ P- 

We can then solve (2.14) to find 

ik{U-{0) - U+{0))z 


M{k) = —A/+(fc) H- M-{k) = |A:| — tanh (|fc|/i+) H-tanh (|fc|/i_) 

P+ P- 


1±{X2) = 


(2.18) 


p±M{k) cosh {h±\k\) 

p^(t7(0)-17+(0)) 


cosh {\k\{x 2 T h±)) 

cosh (|A:|(a ;2 T h±)) 


ikz. 


|A:|(p_ tanh (|A:|h+) + tanh (|A:|/i_)) cosh (/i±|fc|) 
Substituting (2.16) and (2.18) into (2.13) (with the + sign at X 2 = 0), we obtain 

ik{U+{0) - c)z = 3 ^ 2 ( 0 ) - A/+(fc)7+(0) = 3^2(0) - |A:| tanh (|fc|/i+) 7 +( 0 ) 


= 3^2(0)- 


p_(f/_(0) — C7+(0)) tanh (|fc|/i+) 
p_ tanh {\k\h+) + p+ tanh {\k\h-) 


ikz 


which implies 

/p+C/+(0) tanh (|fc|/i_)+ p_{7_(0) tanh (|fc|/i+) 

‘H- p_ta„h(|t|I.O+P^t.„l.(|t|fc-) -d^ = A(0). 

Again we recall that because Y G X(S'), 3 ^ 2 +(0) = 3^2-(0). 

In the next step, we will use the fact that Y is divergence free, along with its boundary behavior, 
to eliminate V and 3^i, obtaining an equation for ( 3 ^ 2 , 2 )- Notice, for instance, that because Y is 
divergence free, 

ikyi + d,,y2 = 0 . 

Thus 3^1 can be determined from 3^2- In light of this observation, the horizontal component of 
(2.12) becomes 

(2.20) - (t/ - c)y!2 + U'y2 - k^U-i + iu' + ikV = 0. 

On the other hand, the vertical component of (2.12) implies 

(2.21) ik{U - c)y2 + ikU-f' + V' =0. 


Using the above two equations, we can eliminate P, obtaining the ODE 

( 2 . 22 ) - 3 ^" + (-^ + k^)y2 + -^ 7 ' = 0 , 

U — c U — c 
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where the prime denotes dx^ and we have used the fact AF = 0. As 1" € X(S'), we have moreover 
that 3^2±(±^) = 0. The behavior of 3^2 on the interface is dictated by (2.19). Recall also that the 
boundary behavior of P is given by (2.15), which, together with (2.20), implies 

{g IpI -aP)z= IpVj 

(2 23) 

= li- Ipiu - c)Xl + IpU'l y2 - e [171 P+7+ + IpU'l'l ), 

as P+1+ = P-1- on S. 

In summary, we have the following result. 

Proposition 2.1. —ike is an eigenvalue of the linearized two phase fluid Euler equation if there 
exist nontrivial linearized solutions of a single Fourier mode of the form (2.16). This is equivalent 
to the existence of a nontrivial solution (z,y 2 -,i,c) to (2.18), (2.19), (2.22), (2.23) and such that 

y2±{±h) = o. 

Remark 2.1. The above calculation is still valid if one or both of h± becomes infinity. For example, 
if = oo, then defining tanh(|A:|/i_) = 1, we get 

, , p+(t/_(0)-t/+(0))el'=l"= 

1 -{X 2 ) = 7777 -; 777771—7":- -tkz, 


y2{0)=ik[ 


\k\{p- tanh {\k\h+) + p+) 
/p+17+(0) + p-U-{0) tanh (|fc|/i+) 


— c]z. 


p_ tanh (|fc|/i+) + p+ / 

Remark 2.2. Note that, up to this point, we have not exploited any small parameters in the problem. 
So, for instance. Proposition 2.1 holds even for a system consisting of two fluids with roughly equal 
density, as one would expect with internal waves moving through a channel. Of course, the heart of 
the analysis to come is in determining (0), and for this we will make strong use of the assumption 
that p+/p- 1. 

Remark 2.3. While U" appears in equation (2.22), it is actually not necessary to assume U'f exists 
in the strong sense. In fact, in the hypotheses of Proposition 2.1, one may replace (2.22) by the 
following equation derived from (2.20) and (2.21): 

( - (c/ - c)y^ + U'y2 - + 7'C/'y = ik{ik{U - c)3^2 + ikUi^. 

This is of course equivalent to (2.22) if U± G , but makes sense even if Uf. has jump discontinuities. 
In the latter case, one expects that will likewise exhibit jump discontinuities at the same locations 
as U'. In particular, Proposition 2.1 is still valid even if [/[_ has jump discontinuities in the bulk of 
the air or water regions, which corresponds to the situation where of U" possesses 5-masses. This 
justifies our consideration of piecewise linear wind profiles in Section 5. 

Remark 2.4. Let us now revisit the question of how our method differs from the formal linearization 
procedure. In (2.7), we demonstrated that daV will not satisfy the correct boundary conditions on 
the interface S unless 

z • {-dxiV*, I)’" - [-c* • ei| dx,z = 0. 

Note that in the special case of shear flow, this simplifies to 

[u* • eil = 0, 

which is precisely the statement that there is no vortex sheet. Hence, for a background flow that 
is continuous over the interface, daV is indeed the right linearized quantity to consider in the sense 
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that it is in X(S'), but the second one allows for a vortex sheet, or a non-laminar flow, this ceases 
to be the case. 

This is not merely a technical point. Observe that the tangential velocity does not affect the 
motion of the fluid interface, and so to truly have a statement about the formation of surface waves, 
one must guarantee that the instability is for the dynamics of the normal velocity. Our splitting 
method is precisely what allows us to do this, and what enables us to see directly what quantities 
must be controlled in order to ensure stability/instability. 

Lastly, we mention that this has an underlying geometric intuition. Consider the Lagrangian 
formulation of the Euler interface problem. Incompressibility is equivalent to the statement that 
the restriction of the Lagrangian flow map to each fluid region is volume preserving. The set of 
such mappings can be viewed as a submanfiold A4 embedded in [Uq] pdx) (cf., [2, 5, 10, 19, 21]). 
Naturally, the linearized problem about a particular Lagrangian flow map uq should then be set 
on the tangent space T„qA 1, and the corresponding variation of the Eulerian velocity should lie in 
uo{TuoM). When one formally linearizes directly in the Eulerian variables, there is no guarantee 
that this will be the case because Ai is not a flat manifold. In effect, by thinking exclusively in 
terms of Eulerian variables, one risks losing a crucial piece of geometric data: the base point of 
the tangent space. Our procedure is carried out in the physical variables, but the splitting is done 
exactly so that Y o uq € TugM. and VT ouq € (TugM)^. 


3. Linear instability and critical layers 

In this section, we consider the physical regime where 

0 < e ^ < 1, C/_ = 0, 

P- 

which means the upper fluid (the air) has much lower density than the lower fluid (the water), 
and that the lower fluid is stationary. Moreover, since we are interested in the problem of wind- 
generation of surface waves, we look for linearized unstable solutions with z ^ Q. Without loss of 
generality, we normalize by taking 

ikz = 1 . 


3.1. Derivation of the dispersion relation. Under the above assumptions, 7 satisfies 

t/+(0) 


7+(0) = - 


(3.1) 


l'+{x2) = |fc|7-s(0) 


|fc|( tanh (|fc|/i+) -I- etanh (|fc|/i_)) 
sinh(|fc|(x2 - h+)) 


, 7 +( 0 ) =-|fc|tanh(|fc|/i+) 7 +( 0 ) 


cosh {\k\h+) 

The boundary conditions for 3 ^ 2 , which solves equation (2.22), take the form 
(3.2) y 2 ±i±h±) = 0, 3 ^ 2 ( 0 ) =-ejfcj tanh(|A:|/i_)7+(0) - c. 

Since U- = 0 in the water, 3 ^ 2 - can be determined explicitly: 

cosh (|A:|(a ;2 -I- /i_)). 


(3.3) 


3^2-(2:2) = 


-3^2(0), 3^^_(0) = |/c|tanh(|/c|/i_)3^2(0). 


cosh {\k\h-) 

Therefore, the unknowns for the linearized systems reduce to (c, 3^2+)- In addition to boundary 
conditions (3.2), 3 ^ 2 - 1 - satisfies 


U" 

-XV + ( + 


U+-C 


k'^)y2+ 


U’l 


U+-C 


7V =0. 
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We first perform a change of variables to transform this equation into a homogeneous one. Let 

_ 3 ^ 2 + + 7 + _ 3 ^ 2 + + 7 + 

3 ^ 2 ( 0 )+ 7 ^( 0 ) “ U+{Q)-C 

Then y satisfies 

(3.4) + -h = 0 on X2 e (0,/i+), y{h+) = {) 

(7+ — c 

along with the normalizing condition 

(3.5) 2/(0) = 1. 

This is simply the classical Rayleigh’s equation that is well-known in the study of the linear in¬ 
stability of shear flows on fixed strips. Notice, however, that the boundary condition for y is 
inhomogeneous due to the interface motion. Returning to (2.23), we see that —ike is an eigenvalue 
of the linearized problem if 

5(1 - e) + — = -e(t/+(0) - c)3^'+(0) - c3^'_(0) + eC/; (0)3^2 (0) 

-efc2c/+(0)7+(0) + eC/;(0)7;(0). 

Substituting (3.1), (3.2), and (3.3) into the above equation we obtain 
erk'^ 

g{l — e) H - = — e([/+(0) — 0)3^2+ (0) + c.[c\k\ tanh(|fc|/i_) — e[/+(0)) 

-I- e|A:|(c|fc| tanh^ {\k\h-) — etanh(|fc|/i_)f7(|_(0) 

- \k\U+{0) - [/+( 0 )tanh(|A:|/i+)) 7 +( 0 ). 

Finally, in terms of y, we have 

2/(1 — e) -f pZ^erk^ = —e(f7+(0) — cYy'{0) -f c^\k\ tanh(|fc|h_) -|- ef7^(0)(?7+(0) — c) 

(3.6) 

ec|A:|f7-|_(0)(l — tanh^ (|/c|h_)) 
tanh (|fc|/i+) -I- etanh (|A:|/i_) 

This is the dispersion relation for the linearized problem with a quiescent ocean and a shear flow 
in the air. 

We summarize the computations above in the following proposition. 


Proposition 3.1. The linearization of the two phase fluid Euler equation at the shear flow v = 
[/+(a;2)xTx[o.?i+]®i eigenvalue —ike if a solution y of (3.4) satisfies (3.5) and (3.6). 


When e = 0, equation (3.6) has a pair of solutions for c: 


(3.7) 


Ck = 


g -I- p_(jkf 
|fc| tanh (|fc|h_) 


or — 


g -I- p_(jkf 
|fc| tanh (|fc|h_) 


This is simply the dispersion relation of the one phase fluid problem (where the air density is taken 
as p+ = 0.) For e in a neighborhood of 0, equation (3.6) is a quadratic polynomial in c with a 
complex parameter ey'flS) and a real parameter e, when we fix others like fc, h±, [/+(0), f7(|_(0), and 
so on. The solution c can be expressed analytically in terms of e and ey'(O) 

c = F(ey'(0),e). 
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Near e = 0 and e?/'(0) = 0, this analytic expression has two branches containing the positive and 
the negative values of Ck in (3.7), respectively. The quadratic formula clearly implies F{a,e) € R 
for small a G R. Therefore, near (e, ey'{0)) = (0, 0), F must take the form 

(3 8) c = F{ey'{0),e) 

= /R(eRe?/'(0),eImy'(0),e) + ie Im?/'(0)//(eRe?/'(0), elmy'(O), e), 
where fa and // are real valued analytic functions satisfying at (0, 0, 0) 

(3.9) «0)=c.. a/„(0) =/,(0) = ft/«(0) = 0. 

Here the formula for //(O) can be obtained via implicit differentiation. More detailed information 
about // and fji can be derived from the quadratic formula if needed. Since we are interested in 
instabilities, we will seek solutions {y,c) of (3.4), (3.5), and (3.6) with Imc > 0 and c near Cfe. 
Clearly the key task is to analyze Imy'(O), whose dependence on c is quite intricate as it involves 
solving the Rayleigh equation with a singularity. 


3.2. Examples. We first present a few examples where the profile is simple enough to do explicit 
calculations. Throughout, h± = oo is assumed for ease of computation. 


Example 3.1 (Kelvin-Helmholtz instability). Suppose that U+ = Uoj he., the wind velocity is 
uniform. Then Rayleigh’s equation (3.4) simplifies to 

(3.10) — y" + k^y = 0 on 12 G (0, 00 ), ?/(0) = 1, y —>■ 0 as 0)2 —>■ 00 . 

This can be solved explicitly. We hnd in particular that 

y'(0) = -|fc|. 

Inserting this into (3.6) yields the following quadratic equation for c: 

g{\ - e) + pZ^crk"^ = e|fc|(17o - c)^ + c^\k\. 

Since the coefficients are all real, instability ensues if and only if there are complex roots of this 
polynomial. Simply evaluating the discriminant reveals that this will be the case if and only if 


(3.11) 




p- - p+ 


+ cr|fc|'’ 


1 


P+ + P- P+ + P- 

Here we have rearranged terms so that the densities ratios are dimensionless. This inequality is 
precisely the Kelvin-Helmholtz instability criterion (cf., e.g., [9]). 


Example 3.2 (Constant shear without a vortex sheet). Consider the situation where 



X2>0 
X2 <0 ’ 


for a fixed p, > 0. This corresponds to a velocity profile which is continuous over the air-water 
interface and has a constant (nonzero) shear in the atmosphere. 

As before, Rayleigh’s equation (3.4) reduces to (3.10), and hence y'(0) = —\k\. On the other 
hand, (3.6) becomes 

y(l — e) + pZ^crk^ = (1 + e)|fc|c^ — e/ic. 

It is completely elementary to show that the above quadratic equation has only real roots when 
0 < e < 1. We conclude that, for any choice of p, the corresponding wave is linearly stable. This is 
in accordance with Miles’s prediction, and our own Theorem 4.1, because U" vanishes identically. 
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Example 3.3 (Constant shear with a vortex sheet). Building on the previous example, let us now 
take ?7+ to be of the form 

U+{X2) ■■= Uq +11X2, 

for some Uq > 0 and fi. The Rayleigh equation for y is trivial to solve explicitly and we find once 
more that y'{0) = —\k\. Thus the dispersion relation becomes 

g{l - e) + —= e{Uo - c)^|fc| + ey{Uo - c) + c^\k\. 

P- 

Evaluating the discriminant, we infer that U is unstable if and only if the following inequality is 
satisfied 

^ TT nr , \ ^ I 5 ^ ~ ^ I 

(l + e)2^°^ ° + + mT+~e + p_(l + e)- 

Rewriting this in terms of p±, we get: 


i2 P+P- 

{p++P-r 


UoiUo + J^)> 


pI 


iP+ + P-)^ 4 


+ 9\k\ 


P- - P+ 
P++ p- 




1 

P++ P-' 


Comparing this to (3.11) reveals that vorticity in the air region — even constant vorticity — can 
be destabilizing in the sense that it may reduce the minimal value of Uq required for the wind- 
generation of water waves. 


3.3. The necessity of critical layers for instability. Formal calculations indicates that there 
exists an unstable eigenvalue —ike with c near Ck provided that Cfc belongs to the range of U+ on 
[0, /i+], see for example [15, 14]. In the following lemma, we prove that this is a necessary condition 
for the existence of instability near Ck- Here, for simplicity, we only consider the case of finite 
atmosphere h+ < oo. 

Lemma 3.1. Suppose U+ € and Ck ^ [/+([0, h+]), then there exists Cfc > 0 such that, if {y,c) 
solve (3.4) along with (3.5) and (3.6) for e G (0,ek), and 

\c-Ck\<\ niin \ck-U+{x 2 )\, 

4 X2G[0,h+] 

then c G K. 


Proof. Let if = ' 0 ( 2 : 2 ) be defined by 

y{x2) =■■ {U+{x2) - c)ip{x2) + 

li+ 

and denote 

S := min jcfe - C/+(x 2 )|. 

X2£[0,h-^] 

Then (3.4) implies 

I -{{U+ - + k^{U+ - cfif + {U'l + P(C/+ - c)) = 0, 

\ 0(0) = ip{h+) = 0. 

Multiplying the above equation by if and integrating on [0, /i+], we obtain 

^ iU+-cfik-^\if'\^ + \if\^)dx2=J^ ( fc2c) {U+-c)ffdx2. 
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Even though c may not be real, |c— Ck\ < i5/4 implies Re (t/+ — c)^ > i5^/2. Taking the real part of 
the above equality, we have 

which implies that 

k-^\y'\L- + \y\L-<C6-\ 
where C is a constant independent of i5, e, and c. 

Let 

(3.12) 

One may compute 

Im c 

(3.13) = W(0) = lmy'(0), W(/i+) = 0. 

|C/+ - c|^ 

Along with the above estimates on y, this implies 

I lmy'(0)| < C5~^\ Imc|. 

From (3.8), we obtain | Imc| < CS~^e\ lmc| and thus Imc = 0. □ 

Remark 3.1. We emphasize that this result, under the assumption on U+, means unstable 
eigenvalues can only bifurcate out of the imaginary axis from ±ikck in the range of U. If U+ ^ , 

as our constructions later show, it is entirely possible that there are unstable eigenvalues lurking 
elsewhere (cf. Section 5). 

The calculation of the Wronskian (3.12) and (3.13) will play an important role in the next section 
where we provide a sufficient condition for instability. 

4. Instability induced by critical layers 

In this section, we present a sufficient condition for linear instability related to critical layers. 
We do this by seeking a solution (y,c) of (3.4) along with (3.5) and (3.6) with |c — Cfe| <C 1 and 
Imc > 0, for e := <C 1. Here Ck is given in (3.7) and we are assuming < oo. Our main 

result is the following. 

Theorem 4.1. Assume U+ € , /i+ < oo, and Cfc S K satisfies 

0 7 ^ {s G [0, h+] I U+{s) = Ck} = {si,..., Sm] C (0, h+), 

and 

U}_{sj) 7 ^ 0, Cfct/"(sj) <0, VI < j < m, and CkU"{sj) < 0 for j = m — 1 or m. 

Fore = P+Ip- ^ 1; there exists a solution {y,c) of (3.4) along with (3.5) and (3.6) with |c—Cfc| = 
0(e) and Imc > 0 with a positive lower bound of order 0{e). 

Remark 4.1. Here Cfc may take either its positive or negative value, whichever satisfies the above 
assumptions. A weaker hypothesis is given later in (4.32) and the leading order form of Imc can 
be found at the end of Subsection 4.4. 

Before we give the rigorous argument, it is interesting to outline the heuristic calculation of 
Imy'(O) which is the key in obtaining the instability due to (3.8). The essence of this calculation 
can be found, for example, in [14]. 

Assume c* G M satisfies 

(4.1) {s G [0,/i+] I C/+(s) = cj = {si,.. .,s™} C (0,/i+), U}_{sj) 7 ^ 0, j = 1,.. .,771 > 1 
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and (y,c) solves (3.4) with |c— c,| <C 1, 0 < |Imc| C 1, and \y\^ reasonably regular. Let 


:= C/+^({Rec}), 

where s' is close to Sj. By integrating (3.13), we first obtain 




where the discrete summation resulted from the limit i5-masses produced by the singularity of the 
integrand at the critical layers. 

From (3.8) we find 



Therefore, instability occurs if the above square bracketed term is negative. We will make existence 


of critical layer induced instability rigorous in this section. The crucial part is the analysis near the 


coefficient singularity of the Rayleigh equation (3.4). 

Rayleigh equation In the process of proving Theorem 4.1, we obtain the following proposition on 
the convergence of solutions to the Rayleigh equation (3.4) as the parameter c approaches a limit 
in R. This result can be useful by itself in the study of instability of shear flows. 

Proposition 4.1. Suppose c* G R is a regular value of G C\ I > A, on [0,/i+] and c* 
{f7+(0), f/+(/i+)}. We use the notation as in (4.1). For c = + ici G C\M sufficiently close to 

c*, let y be the solution of (3.4) along with y'{h^) = 1. We have 

(1) There exists a unique solution y* of 



that, at any s G f7_|_^({c*}), exhibits the behavior 



|yl(s + X 2 )\ = 0(1 log 1 x 211 ) as X 2 -)■ 0. 


Moreover y„ is ^ in c* for X 2 ^ 0^^({c*}). 
(3) For any 5 > 0 and a G (0,1), 


m 


\yix 2 ) - y*{x 2 )\ = 0(|c/|“), for all X 2 G [0,/i+] \ |J(sj - S,Sj +5). 


The above estimates are uniform in k for k on any compact subset of . 
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Outline of the proof of Theorem 4.1. While the above formal argument provides a useful 
insight, it is far from straightforward how to turn it into a rigorous proof of Theorem 4.1. Among 
the issues, for example, are that one needs some control on y and y' for |Imc| 1 , including 
some positive lower bound on |yp to ensure the instability. This can be potentially achieved by 
identifying the limit of y as Imc —^ 0 along with convergence estimates, but this is very nontrivial 
due to the creation of i5-masses near critical layers in the limiting process of the singular equation 
(3.4). As the rigorous proof of Theorem 4.1 presented in this section is rather technical, here we give 
a brief outline of the analysis near one singularity sq G (0, /i+) where |17+(sq) — Rec| + | Imc| <C 1. 

Step 1. As Imc —>■ 0, and y' do not remain bounded uniformly. Our first step is thus to un¬ 
derstand the behavior of y{x 2 ) for X 2 near sq when |Imc| <C 1. We change variables to 
T := U{x 2 ) — Rec, which is more convenient for the local analysis due to its appearance in 
the denominator. Meanwhile the unknowns are transformed to the more geometric quanti¬ 
ties (ui, U 2 , U 3 , W), where ui and U 3 represent the squares of the norms of y and j/', and U 2 
and W the dot and cross products of y and y' . It turns out that ui remains uniformly Holder 
continuous as Imc —>■ 0, while W develops a jump discontinuity at r = 0, U 2 a logarithmic 
singularity in r, and M 3 a singularity of the order of the square of logarithm. More care¬ 
ful analysis reveals that in the leading orders the singularities are symmetric in r near t = 0. 

Step 2. A priori estimates in Step 1 motivate us to make the right guess for the limit equation 
(4.12) along with conditions (4.13) at the singular point t = 0. Coefficients and solutions 
of this limit system still possess singularities. To better understand the limit problem, we 
apply another linear transformation B{t) to the unknowns which depends on r smoothly 
in neighborhood of r = 0. The resulting system has very simple variable coefficients and 
can be solved explicitly. 

Step 3. After carefully separating the singular parts of the solutions, and with the help of the above 
linear transformation B{t) applied to solutions y of (3.4) for |Imc| -C 1, we will complete 
the proof of the Proposition 4.1 and obtain very good error estimates near one critical layer. 
Finally, the proof of Theorem 4.1 is carried out by combining these estimates near all of 
the critical layers. 

4.1. Preliminary estimates near a singularity of Rayleigh’s equation. Assume U+ G ^^([0, /i+]) 
and suppose there is an sq G (0, h+) and <5 G (0,1) with 

(4-3) t 7 +(so) 7 ^ 0 , g (1^2) on [so-<5,so + <5]. 

In this section, we analyze the solutions y of (3.6) on the interval [sq — sq + for c very close to 
[/+(so) with c/ := Imc 7 ^ 0. For such c, there exists s such that 

Rec =: Cfl = 17+(s), and |s — so|, |c/| <C <5. 

In other words, s is the altitude at which the wind profile takes the value Re c. 

Since (3.4) is invariant under multiplication by a complex scalar, we make the following trans¬ 
formation which better reflects the rotational invariance on the complex plane: let 

1 % 

Ml := u2-.= -{y'y + y'y), U3 := |2/T> W :=--(yy'- yy')- 


(4.4) 
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Geometrically, ui and its are the squares of the norms of y and y', while U 2 and W are the dot and 
cross products of y and y'. Obviously, they satisfy 

(4.5) {U 2 f + - Wilts = 0 


and determine the solution y and y' up to a rotation in the complex plane. One of the advantages of 
these new variables is that identity (3.13), which relates W to |i/p, is easier to handle and directly 
links to the generation of the instability. 

One may compute that (iti, it 2 , ws, W) satisfy 


(4.6) 


u'l = 2 lt 2 

Ito = I + 


= 2\e 


w = 


K iU+ - cr) 

{u+ - cnf + c] 

^ UI {U+ - ch) ' 
{U+ - Cnf + cj^ 
ciUl 


Ml + Ms 


M2 


2ciU'l 


{U+ - Cr) + c] 


-w 


{U+ - Cr) + C 


■Ml 


where ?7+, U" are evaluated at X 2 - It is straightforward to show that (4.5) is satisfied by 
solutions of (4.6). 

To handle the singularity U+ — cr = U+ — t/+(s) at s = 0 for small c/, we work with a new 
independent variable 

T = r(a; 2 , s) := U+{x 2 ) - cr = U+{x 2 ) - U+{s). 

This is a valid change of coordinate on [sq — 5, sq + <5] depending on the parameter s. Let 
<5 i.2 = <^ 1 , 2 ( 5 ) > 0 be defined by 


—^i(s) := min{t/+(so ± i5) — f7+(s)}, ^ 2 ( 5 ) := niax{C/+(so ± <5) — U+{s)}. 


Note that this implies that 81^2 = 0{8) due to (4.3) and the fact |s — So| 8 . For r £ [—(5i, ^ 2 ], we 
have 


(4.7) 


where 


MIt = 71-M2 

01 

w, = T_g2^„. 

Ui r2 + ci 


TU 2 


+ Cj 

TU 2 
+ cj 


Ml + 


M2 + 


Ui 


2ciU2 


(r2 + cj)Ui 


W 


Ui{s,t) := f/+(a; 2 (r,s)), U2{t,s) := G"(x 2 (r, s)). 

Clearly Ui ^2 and (5i,2(s) are also smooth in s. 

We first obtain some uniform estimates on solutions to the transformed system (4.7). 


Lemma 4.1. For any a £ (0,1) and Cq > 1, there exist 8 £ (0,1) satisfying (4.3) and C, eo > 0 
depending only on a, k, Cq, |t/(|_(so)|~^; Mwd |t/+|c3([o,/t+])7 such that for all 

|s - Sol < ,0 < \ci\ < Co and ti _2 S [- 81 , 82 ], with \t 2 \ < CqItiI, 
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solutions to (4.7) satisfy on the interval [Ti,r 2 ] (or [T 2 ,ri] if T 2 < ti) 
|wi(t)|, |>V(r)| < C\u{ti)\^ 


|^^ 2 (T)| < C{1 + I log 


-I- 
C? + 


|M(ri)U, |w3(r)| < C'(l + (log^^|-^)2)|u(ri)|. 


\ui{t') - Ui{t”)\ < -t"\+ {cj + {T'"f)^\T' -t" 

where t'" > ^max{|T'|, |r"|} and 

W{t)\w ■= (mi(e)^ + M2(t)^ + usirf + >V(t)^) ^. 

Moreover, for t < min{|Ti|, |t 2 |}, 

|m2(t) - U2[-t)\ < Cr^^““T“|u(ri)|u,, |m3(t) - U3{-t)\ < C\u{ti)\,, 


-/ 

1-0 


i(n)U, 


where 

r. = (c? + r2)i 

In the above estimates, an almost logarithmic singularity appears at r = 0, but disappears as r 
evolves past 0 at roughly the same rate at which it appeared. 


Proof. For notational simplicity, we mainly consider the case ti > 0 and T 2 € [—ti, ti], the argument 
for the other cases being similar (or easier) due to the absolute value outside the logarithm in the 
above inequalities. See the comments at the end of the proof. 

First, we establish an inequality that we will make repeated use of later: 


(4.8) 



log 


cj + (T"'f 
cj + 


dT<C(cf + (T'"f)^\T'-T"r, 


where |r'|, |r"| < Coir'"!, 


and the constant C depends only on a G (0,1), Co > 1, and m > 0 but is independent of c/, t',t", 
and t'" . The above estimate is based on the observation 


cf + ^ 1 

:? + r2 -Cl 


CJ + T 

Applying Holder’s inequality 


, cj + ir"r 


< C 


cj + {T"'f 


cj + 


log 


cj + {T'"f 


cj + r"’ 


a r" 

T~^dT^ 

< C{cj + (T"')")^|r"'|^|r' - r'T 


1 —Q! 


and thus (4.8) follows. 

For a constant M > 1, let 

To := sup{r' G [r 2 ,ri] | |u 2 (r)| < M{l + \ log |)|u(ri)|^,, for all r G [r',Ti]}. 

Clearly tq < ti. We will show tq = T 2 for appropriately chosen M and 5. In the rest of the proof, we 
will use C to denote a generic constant depending only on Co, ct, k, |C/)|_(so)|~^i and |[/+|c3([o,?i+])- 
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For tq < t' < t" < Ti, let t'" be given with t'" > ^max{|r'|, |r"|}. Then from (4.7) and (4.8), 
we have 


\ui{t')-ui{t'')\<CM\u{ti)\^ f (1 + I log ^ \)dT 

Jt' Cj + t 


< CM\u{Ti)\yj(^{l + I log 


cf + tI 


W - t"\ + 


< CM 


(|log 


- t"| + (c? + (r"')')^ 


c] + (r "0 


CJ + T-^ 


This inequality indicates that ui is Holder continuous. It also shows that |ui| < 2 |m(ti)|u, on [tqjTi] 
if M and S are chosen to satisfy C6M < 1. 

Integrating Wt on [r, ti] with r S [To,ri] and using the estimate of ui, we have 


|W(t) - W(ti)| = I j ^ (r0^ + ~ 

if CSM < 1. 

We will proceed in two steps to estimate ^2,3. First, for t £ [max{0, tq}, n], (4.7) and the above 
bounds on W imply 


|M3(r) - M3 (ti)| < C\u{ti)\^ 


M{1 + 




r) (1 + I log 


l) + 


ci 


ci + Tf 


c2 + (^/)2D (^,)2 + 


dr'. 


Some terms in the integrand above can be directly integrated, yielding terms like tan log, and 
log^. Ultimately, we hnd 


|u 3 (t) - U3{ti)\ <c(^1 + M{ 


n In 


+ I log -TT^I + (^°g *^2 t ^2 )^)) 

ci + M ci + M / 


c? + T? 


Along with the estimates on ui_ 3 , the U 2 equation in (4.7) implies, for r £ [max{0, tq}, ti], 


|m2(t) - U2(ri)| 


(1 


cj + (r')2 


7^)|u('ri)U + \Mt')\ 


dr' 


< Cdlog^^l +Mri-“|ri -rr)|«(ri)|, 
Ct ~r r 


C? + T? 


< C{1 + |log i ^J )|m(ti)|^, 

Cj ~\~ T 


if CMS < 1. From the above inequalities, we infer that the value M as in the definition of tq can 
not be achieved on [max{0, tq}, ri] if M is reasonably large and 6 is chosen such that CSM < 1. 
Therefore, either (i) tq = T 2 in the case T 2 £ [0,ri] which completes the proof of the claim tq = T 2 , 
or else (ii) tq < 0 in the case T 2 < 0 . 

Let us consider the latter: suppose that T 2 < 0. Unfortunately, one can not simply extend the 
above argument for r < 0 since the desired logarithmic upper bounds on 1 x 2,3 start to decrease as 
T decreases past 0. Instead, we study the quantities 


M2,3('r) := W2,3('r) -'«2,3(-'r), T > 0, 
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which satisfy 


(4.9) 


U2t = 


^3t — 


1^.1 Is 

+ (tt- + 777—w )“3 + 


,U 2 Ui U2U1 


Ui{-t) 'C/i C/i(-r) 

k'^ui k^ui, . 

-(-r) 


+ cf ^ Ui 


Ui 


(-t)) 


Ui Ui 
2t /U 2 , 


t2 + C? 


'El-El 

'Ui Ui 


{-t))uE) +2fc2(^ + m(_^)) 


, 2c/ ,,C/2W, , ,C/2W,, ,, 


Here, all functions are evaluated at r unless stated otherwise. Recall that the assumed upper bound 
on U 2 , and the estimates derived for ui and W, are valid on [ro,Ti] and tq < 0. Using these, with 
t'" = |ri| for ui, we obtain on [ 0 , min{— tq, ri}] 


(4.10) 


\U 2 \ j (|M 3 (r')| + |m 3 (e')| + [MtI “(r')“ ^ + l)|u(ri)|u,) dr' 

Jo 


< 


C ^MtI “t“|m(ti)|u, + j |u3(t')| dr'^ . 


Similarly, the M 3 equation and the assumed upper bound on M 2 imply that, for r G [0, min{— tq, ti}]. 


< c 


0 Cj + {T') 


yw|M2(T')|dT'+ (l + Mri “r“)|M(ri)|i 


Substituting the estimate (4.10) for M 2 into the one for M 3 above, taking CMS < 1, and then 
integrating by parts, leads to the following inequality: 

l^3l <C'(|M(r3)U+ rMri-»(TO“-V(ri)U+ , f \u^{T")\dr"dr' 

^ Jo + V ) Jo 

<c(\u{ti)\^+ f |log^^^-i-^||M3(r')MT'). 

^ do + (T ) / 

Therefore, for r G [0, min{— tq, ti}], 

|M 3 |cO([ 0 .r]) < C (^|m(ti)|„ + | log | dr'^ 1^3 | ([O.r]) 

< C(|M(ri)|™ + (cj +T^)^r“|M3|cO([0,r])), 

which, along with the fact that M 3 ( 0 ) = 0 , implies 


|w3|co([0,r]) < C'|M(ri)|^ 

for T G [0, min{—T q, Ti}]. From this we may conclude that 

|m2(t)| < t G [ 0 ,mm{-To,Ti}]. 

Now, choosing M and 6 such that M > C and CMS < 1, the above inequality and the previous 
estimate of U 2 on [ 0 ,ti], show that we must have tq < max{— Ti,r2}. Thus, all of the inequalities 
in statement of the lemma hold on [max{T 2 , — t/}, t/], which completes the proof in the case T 2 G 

[-Tl,Ti]. 
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We will conclude the proof of the lemma by discussing several other cases of ri_ 2 - First note that 
the case ri < 0 can be treated in the exactly same manner as the case ri > 0 because one may 
consider u{—t) and all the estimates go through. The above arguments complete the proof in the 
case Ti > 0 and T 2 € [—In fact, if T 2 € [0,ri], the estimates on { 12,3 are superfluous. The 
proof for the case 0 < ti < T 2 < CqTi is exactly the same (again, there is no need to consider { 12 , 3 }■ 
Finally, if T 2 € [—CoTi,—ri] we can combine the estimates on [t 2 ,—ti] and those on [—Ti,ri] to 
obtain the desired inequalities. □ 


The above lemma provides some preliminary control of solutions to (4.7). We refine them in the 
next corollary. 

Corollary 4.1. For any ri _2 G [—(5i,^2] with |r 2 | < |ti|, the following inequality holds 

|ui(ri) - Ui(t2)| < C(|C/| + log^(|c/| + |ti|)|ti - T2\°'\u{d2)\w 

Moreover, for 0 < r < min{i5i,52}, we have 

\u 2 {t) - U 2 {-t)\ < Ci\ci\ + |T|)i-“r“log2(|cz| + |r|)|u(<52)U 

|W(r) - W(-r) - tan-i -| + |W(r) + W(-r) - 2W(0)| 

Fi(0) Cl 


< C\ci\ log^ \ci\ l0g(l + r-r)|u((52)|u,, 

c/ 


and 


|M3('r)-W3(-T)-^^(W(-T)+W(r))tan ^ < Ct“|m(,52)U. 

Here C is independent of ci and r. Also, wherever |m((52)|u, occurs above, it can be replaced by 
\u{-Si)\w 

Proof. Taking t'" = |ti| and appealing to Lemma 4.1, we can see that ui is Holder continuous on 
the interval [T’i,r 2 ] with exponent a and Holder constant 

C{cj + rf)^|u(Ti)|u, < C(|c7| + log^(|c/| + In|)| m(i52)U 

where we used the estimates on |m(t)|i„ in terms of |u(( 52 )|m. This proves the improved Holder 
estimate of ui in the corollary. The estimate on U 2 (t) — U 2 {—t) follows similarly. From (4.7), we 
have 

5r(W(T) - W(-t)) = 72 ^((^^)W + 

Thus, the Holder continuity of ui implies 




Cl 


c/ + r 


■log^dc/l + r) |m((52)U. 


For c/ ^ 0, W is smooth and so we are permitted to integrate this identity yielding 


W(r) - W(-t) - 2^^^(0)tan ^ ( —) < C|c/||m((52)|i 
Ui Cl 


log^ (|cj| +r') 
|c/| + r' 


<C|c/||M(^ 2 )Ulog^ |c/| 


1 


-dr'< C|c7|log |c/|log(l + —)|u(( 52)U, 
/o |c/|+t' |c7| 
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and thus we have proved the desired estimate for W(r) — W{—t). Again from (4.7), we have 

(4.11) 


laAW(T) + W(-r))| = |-A^((i^)(,) - (i^)(-r))| 


J 

2 / 


Ui 


/ ^ c/ log ( C/ + r , 

< C -M <l2)k 

T + C/ 


and the desired estimate on >V(r) + W(—r) follows similarly. 

Finally, to control U 3 (t) — it 3 (—r), we use (4.9) and (4.11), along with Lemma 4.1 (first with 
Ti = T and then ti = 62 ), and the above estimate on W, to show 


\dr(u 3 {T) - U 3 {-t) - ^^(W(-t) + W(r))tan ^ (^))| 

< C(log^(|c/| +t)\u{S2)\w + 2 ] 2 I^ 2 (t) - U2(-r)| + |m(t)|u,) 

Cj “T T 

< C(log^(|c/| +r)|it(d 2 )U + I“(t)U) < C\o^{\ci\ +T)|u(d 2 )U. 


Integrating this inequality using (4.8) yields the desired estimate for U 3 . 


□ 


4.2. Heuristic limit equation near a singularity of Rayleigh’s equation. Consider a family 
of solutions (u = (ui, U 2 , ms)^, W) parametrized by c/ ^ 0. As c/ —>■ 0±, formally, we expect (u, W) 
to converge to a solution {V = (Vi, V 2 , Vs)"^, W*) of 


(4.12) 


= 


V Ui 


—V 2 


t( 


3 


=: iA(r)C, 


W4.r =0 


on the interval r G [—di, (52]\{0}. At r = 0, Lemma 4.1 and Corollary 4.1 hint that 

• Ui has a continuous limit; 

• lim^^o limcj^o(M2(T) - U2{-t)) = 0; 

• W converges to a piecewise constant function with a jump (Q) at t = 0; and 

• limr^o limcj^o(M 3 (T) - U 3 (-t)) exists. 

Therefore, we impose the following conditions at the singularity r = 0: 

(^i. 2 (r) - Vi, 2 (-r)) = 0, W 40 +) - W 40 -) = 

(4.13) ^ 

(^ 3 (t) - Vsi-r)) = sgn{ci)^^^^{yV40+) + W^O-)). 

Even though c/ is already taken as 0 here, we still track its sign as the direction from which the 
real axis is approached does make a significant difference. 

We begin by analyzing the limiting system which will help us justify the convergence and obtain 
the necessary error estimates. Since V is decoupled from W in the limiting system except at r = 0, 
most of the work can be first carried out for V. 

As a preliminary step, in the following lemma, we perform a coordinate change that makes the 
1/r singularity in (4.12) more tractable. 
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Lemma 4.2. Assume U+ G , I > 3, then there exists B[t) G M 3 x 3 (K) (Msxsi^) the space of 
3x3 real valued matrices) which is in t G [— 61 , 62 ] such that B{0) = I and for any solution 

V of (4.12) for T ^ 0, V ■= B(t)~^V satisfies 


(4.14) 

where 


K = 


t/ 2 ( 0 ) 

tC/i(0) 


AqV = t-^A{Q)V, t^O, 


/O 0 0> 

Ao = 1 0 0 

\0 2 Oy 

Moreover, there exists j3 G (0,1) such that B is dependence on the parameters k and 

s for s is close to sq. 

It is worth pointing out that B is smooth at r = 0 even though the equation is singular there. 

Proof. The proof of the lemma is an application of the invariant manifold theory in dynamical 
systems. Substituting V = B{t)V into (4.12), we obtain 

Vr=T-^B{T)-^{A{T)B{T)-TB'{T))V. 

Comparing this to (4.14), we see that it suffices to find B G so that B{Q) = I and 

(4.15) B'{t) = T~^ [^A{t)B(t) — i?(T)yl(0)). 

To remove the singularity r“^, we treat r as a new dimension in the phase space and consider 
a nonlinear augmented system in a new independent variable 9\ 


(4.16) 


A{t)B - BA{0)\ 


=:GiB,T). 


Clearly solutions of (4.16) with r 7 ^ 0 lead to solutions of (4.15). 

Notice that (/,0) is a steady state of the augmented system (4.16) and 


DG{I,0) = 


A AriO) 
0 1 


where the operator A is defined as AB = A{0)B — BA{0). 


One may verify through direct computations that 

^9AB = ^eAiO)Be-SAiO)^ 

Since 

/ 1 


(4.17) 


^ gM( 0 ) ^ 


1 / 2 ( 0 ) 

C/i(0) 

u2{or 


e 


0 

1 

21 / 2 ( 0 ) 


0 \ 

0 


0 1 


V 1/1 (0)2 c/i(o) 

we infer that the solution map steady state (/, 0) has the 9-dimensional invariant 

subspaces {r = 0} C and has only algebraic growth in 9 there. In the transversal 

direction, clearly 1 is an eigenvalue of DG{I, 0). Motivated by the above formula of one may 
verify through a routine calculation that (Bi, 1 )^ is a eigenvector of the eigenvalue 1 where 


Bi = 


n + oo f+oo 

/ e-'^+^^Ar{0)d9 = / 

7o Jo 
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The convergence of this integral is guaranteed by the fact that has only algebraic growth. 

Therefore, from the standard invariant manifold theory in dynamical systems [ 8 , 7], there exists 
a one-dimensional (7^“^ locally invariant unstable manifold of (/, 0 ) which is tangent to {Bi, 1 ) at 
(/,0). The regularity of the unstable manifold follows from the observation that G G C^~^. The 
smoothness of the unstable manifolds with respect to external parameters — so and k here — 
is like the smoothness of the unstable fibers with respect to base points in the center manifold. 
Therefore it is for some /3 € (0,1) [7, 3]. The unstable manifold can be parametrized by r 

as B = B{t), and will satisfy B{0) = I, B G . The invariance of the one-dimensional unstable 
manifold implies that B{t) is a solution of (4.15). Since (4.15) is a linear ODE system, the solution 
B(t) can be extend to r € [— ( 5 i,( 52 ]- In addition, the tangency of the unstable manifold at {1,0) to 
(ill, 1 ) implies 

f + oo 

(4.18) B'{0)= e-^e^^^°'>Ar{0)e-^^^°Ur 

Jo 

which can be integrated explicitly. This complete the proof. □ 

As a corollary, we obtain explicit forms of V (r) = B{t)~^V{t) for any solution (V, W*) of (4.12) 
and (4.13). 

Corollary 4.2. Let (E, W*) be a solution of the limiting system (4.12) satisfying the conditions at 
the singular point in (4.13), then 

(1) V is uniquely determined by constants ( 01 , 02 , 03 ) along with either >V*(0—) or W*(0-|-) 
such that V{t) = B{t)~^V{t) is given by 

Ei(r)=oi, E 2 (r) = oi^^log|r|-^ 02 , 

^3(t) = -f oaiog 

where the ± is chosen so that ±t > 0, and 

03 =sgn(c 7 )^^|^(W 40 -f)-f W,( 0 -)); 
moreover, Vi(0) = Vi(0) = oi; 

(2) there exists a constant C > 0 depending only on k, |t/+|c 3 ([o,h +])7 |C^-i-('So)|~^; ood S such 
that 

i(|E(-di)| + |W.(-<5i)|) < \V{52)\ + |W.(<52)| < C'(|E(-di)| + |W.(-<5i)|) 

for any solution of (4.12)-(4.13). 

Proof. From Lemma 4.2, there exist (oi±, 02 ±, 03 ±), such that 

V{r) = e('°8l"l)^(°)(oi±,02±,03±)^, ±t > 0. 

Since B{0) = I and B{t) is in r, we have for 0 < r -C 1, 

V{t) — V (— t) = B{t)V (t) — B{—t)V (— t) = V{t) — V (—r) + 0(|t| log^ |r|) 
which implies 

lim (V(t) — V(—t)) = lim (V(t) — V(—t 

T-^0+ ^ n r->-0+ ^ ' 
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From the boundary condition (4.13) and the explicit form (4.17) of it is easy to see that 

ai_2+ = ai,2- and the jump condition on a3± = 03 ± 03 holds. For Vi(0), since B{t) = / + 0(|t|), 
we have 

1 / 1 ( 0 ) = lim Vi{t) = lim(l/i(T) + 0(|r| log^ |t|) = oi. 

T^O r—)-0 

Part (2) is a simple consequence of Lemma 4.1 and part (1) of this corollary. In fact, one can 
first identify parameters 01,2 as well as (03 — a3,yV*(0—)) by the values of (1/, W*) and B at —5i. 
Consequently, condition (4.13) allows us compute the values of >V*(0+) — W*(0—) and thus Og and 
og using Vi(0) = Vi(0) = oi. These parameters completely determine the solution on the interval 
[—<51,^2] • Therefore, the estimates on V and W* follow from the explicit forms of V and W* and 
the smoothness of B. □ 

Before we finish the analysis on the solutions to the limit system (4.12)-(4.13), we consider the 
conservation law (4.5). Observe that, at this stage, we cannot yet conclude + W* — I/1V3 = 0 
because (P, W*) is only assumed to be a local solution for r € [—i5i, ^2]- 

Lemma 4 . 3 . Let (P, W*) he a solution of (4.12)-(4.13) for t G [—<^1,(52] • 

(1) Conservation law (4.5) holds in the sense that 

(P2)^ + W* — Pi4^3 = const on [—i5i,52]. 

(2) Suppose (P2)^ + — P1P3 = 0 and Pi(0) = 0, then (P, W,) is smooth on [—5i,52] o,nd 

takes the form 

P(r) = B(r)(0, 0, og)^, W,(0-) = W,(0+) = 0. 

Proof. From (4.12), it is easy to verify through direct computation that dr{V 2 + W* — P1P3) = 0 
for r ^ 0, and thus we only need to show 

{Vi + w! - PiP3)(o-) = {Vi + _ 141/3) (0+). 

Again since B{t) = I + 0(|t|) implies P(t) = P(t) + 0(|t| log^ |t|), it suffices to prove the above 
identity for (P,W*), which is straightforward from the explicit form of P given in part (1) of 
Corollary 4.2. 

Suppose Pi(0) = 0, then Corollary 4.2 implies oi = 0 and thus 

P(t) = (0, 02, log |r| + 03 ± a'sf. 

The smoothness of B and the fact that B{0) = I then give 

P(t) = (0,a2,^^^^H^log|r|+a3±a3)^ + O(|rlog|T||), for |r| < 1. 

Substituting this into the assumption Vi + = P1V3, we obtain 

al + Wl+0{\T\og\T\\) = 0(|r|log^ |t|). 

From this we conclude that 02 = yV*(0—) = yV*(0+) = 0, and hence part (2) follows. □ 
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4.3. Convergence estimates near a singularity of Rayleigh’s equation. Operating under the 
assumptions in (4.3), we now derive an error estimate comparing solutions of (4.7) and the limiting 
system (4.12)-(4.13) for 0 < |c/| <C 1. The following lemma gives some preliminary bounds. 

Lemma 4.4. Given any solution {u,yV) of (4.7) on [—< 5 i, 52 ] such that |u(— i5i)|u, < 1, let (V, W..,') 
he the unique solution of (4.12)-(4.13) so that 

W,(0-) + W,(0+) = 2W(0), C(-<5i) = u{-6i). 

Then we have 

\ui{t) - Ci(t)| < C|c 7 |log^ |c/|, T e [-(5 i,(52 ], 

and 

|m( 52) -m)| < CIc/r, |W( 52 ) - w,(0+)| + |W(-^i) - W,(0-)| < C|C 7 | |log |c 7||3 . 

Or, alternatively, we may assume |u(( 52 )|io < 1 and the initial condition V_{52) = which will 

give the same inequalities except with S 2 replaced by —(5i. 

Note that in general u — V_ does satisfy the error estimate on the whole interval [—( 5 i,( 52 ] as 
Ml —y_i does. This is natural since M 3 and >V are smooth for 0 < |c/| while and are singular 
at r = 0 . 


Proof. Since the linear transformation through B{t) simplifies the equation for V, which is basically 
the principle part of the equation satisfied by (m, W), we apply the same transformation to the latter 
by letting u = B{t)u. One may compute 


(4.19) 



U2{0) 

tUi{0) 


AqU - 


C/C/ 2 M 1 


(r2 + cj)Ui 


C/D 2 


t{t'^ + cj)Ui 


B ^AqBu + B 


-1 



2 c/[/ 2 W 
(r2 + cf)t7i 


T 


where B and f/i _2 are evaluated at r unless otherwise specified. 

Ideally, one might try to control the error m — C = B{t)~^{u — V_) and W — W.... However, M 3 
and W are continuous at r = 0 for c/ 0 , while Hg and are not, and this discrepancy makes 
it inconvenient to estimate m — H and W — W.,, directly. Instead, we consider the modified error 
functions {Z = {Zi, Z 2 , Z^Y','W) defined for r > 0 by 


2'i, 2 := '«l,2 - Zl,2 - {y‘l,2{-T) - Zi,2(-t)) 

= Ul,2 - iil,2(-T) 

2t/2(0)W,(0+) . _1 ' 


^3 := U 3 — V^ — 


Ui{0) 

- ( M3(-'r) - H3(-r) - 

2 C/ 2 ( 0 ) 


(tan i( —) --sgn(c/)) 


2t/2(0)lV,(0-) 

[/i( 0 ) 


(-tan 1 (—) +-sgn(c/)) 

Cl Z 


= Mg - U 3 (-t) - 


(W.(0D+ffi.(0-))..n-R-) 


w := W - W,( 0 +) + W(-r) - W,( 0 -) 


= W + W(-t)-2W(0), 


(4.20) 
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where all functions are evaluated at t unless otherwise stated. Here we have used the property that 
Vjir) — r) = 0 for j = 1, 2, and, for j = 3, is a specified constant given by Corollary 4.2. For 
any c/ 0, {Z, W) is smooth for r > 0 and 


|Z(0)| = |W(0)| = 0. 


From Corollary 4.1, we have 


|W(t)| < C\ci\ log^ \ci\ log(l + -p-r). 

|c/| 

The next step is to estimate Z{t). Towards that end, we compute that 

Zr{T) = Ur{T)+Ur{-T) - (o, 0, (>V J0+) + >V„ (0-))) . 

Recalling (4.19), this becomes 


-(^H ^{Q,ui,2u2f){-T)) +^^(^{B 1(0,1 


^C/i 
- ( 0,0 






um 


um 


Since Ao(0,0,1)^ = 0, we see that, for t >0, Z solves 


+ c/,(0KA+*»;)*°’ 


where (pi ^2 are given by 


Mr) ■■= ^2 ^^2 ^(0,0, ^^)'^)(r) + (H ^(0,0, (-r) 


-(0,0,^(W(r)+W(-T)))^) 
= 0(|c/|r(r^ + cj)"i). 


and 


Mr) ■■= - ^2^ {^i^ ^(0,Mi,2M2)^)(r) - i( 0 , ui, 2^2)'^) (-r)) 

= c? 0 (r“-i(r + |c,|)-(i+“) log'dc/l + r)) = c] 0 {t^-\t + |c/|)-(i+“) log^ \ci\). 

Here we used the smoothness of B and Corollary 4.1. 

From the equation for Zt and the fact that e'’^“(0, 0,1)^ = (0, 0,1)^ due to (4.17), we have 


+ e^(‘°®^)^“((/)i(T') + </> 2 (T'))l dr'. 
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Employing the estimates of (j)i ^2 and W derived above, as well as the explicit expression (4.17) for 
we then Hnd that for r > 0 

T 7'^ rc? log^ |c/| log(l + T^) 

|Z(r)|<C(l + |log-n|Z(ro)|+C / ^ ^ ^2 

, cflog^lc/l 


+ ( 1 +iiog—n( 


C/ T 


dr'. 


'(|c/|+t')2 (t'Y “(|c/| + 

When c/ 0, Z is smooth for r > 0 and thus Z(to) = 0(|to|) as tq 0+. This implies 


\Zir)\ <C r 
Jo 

< C\ci\ [ 
Jo 


cj log^ |c/|log(l + |f^) 

{\ci\ + r'Y 


(1+iiog^n( 


|c/ r 


cj log^ |c/| 


\cj\+t'Y (r0i-“(|c7|+r')i+“' 


dr' 


TvT r' + log^ |c/| log(l + t') 


log^ |c/| 


Y + r'Y 


+ T2(|c/|r') 2( 


log^lwl 


'(l + r')2 {t'Y “(l + r')i+“ 

Integrating the above inequality gives 

\Z{t)\ < C\ci\{\og^ |c/| +log(l + 1 ^)) + Crtlc/p-t log^ |c/| < Clc/I", 
for 0 < r < (5o := min{(5i, <52}. Recalling the definition of Z, and using the fact that 

tan-i( —)-sgn(c/)^ =0(|cj|), 
c/ 2 

we can turn the above inequality into an estimate for u and V_: 

|m(^o) - m)l < \u{-So) - r(-Jo)l + c\ciY. 

Since the system (4.7) is a regular perturbation of (4.12) on [—i5i, 52]\[—i^o, ^o], with the difference in 
their coefficients being of order 0 (|c 7 |), and noting that we have already established that V_{—5i) = 
u{—5i) and B{t) is bounded, it is simple to show that the above inequality implies the desired 
estimate on u{52) — Y_{S 2 )- 

Next, we refine the bound on ui. Observe that the first row of the matrix Aq vanishes. Therefore, 
for T G [—(5i, ^ 2 ], equation (4.19), the property B{0) = I, and Lemma 4.1 together imply 

Ccj log^dc/l + |t|) 


\uir\ < 


Ccj 
cj + t'‘ 


|m(t)| < 


cj + 


L-j ~r / 

Arguing as we did before when we integrated Z^ in the previous paragraph, we obtain 

\ui{t) - £ti(-5i)| < C|c/| log^ |c/|. 

From Corollary 4.2, V_i is a constant and thus C(—^i) = {i(—(5i). These facts, together with the 
boundedness of B and the above inequality lead to the stated inequality for ui — V)j. 

Finally, we complete the proof of the lemma by deriving the estimate on W(52). In fact, the 
bounds on W in Corollary 4.1 give 

- ^7^(0) tan-i ( —) - 2W(0)| < C\ci\ log^ |c/| log(l + -^). 

Ui Cl c/ 


(4.21) 


|2W(r) 
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Moreover, the definition of W,, and the jump condition on at r = 0 given in (4.13) show that 

|W(t) - W,(0+) - ^(0)(Mi(0)tan-^ (—) - sgn(c7)|yi(0))| < C\ci\\og^ jc/l log(l + y^). 

Dl Cl 2 ' \ci\ 

Thus the estimate on ui — V_i implies the desired estimate on yV{S 2 ) — VV,f(0+). The analogous 
inequality for >V(—(5i) — VV,f(0—) is obtained similarly. □ 

Notice that VV,, in Lemma 4.4 does not have the same initial value as W at t = —(5i. This can 
be fixed easily by using Corollary 4.2, as we show in the next corollary. 

Corollary 4.3. Given any solution (u, W) of (4.7) on [— i5i,(52] such that |u(—5i)|uj < 1, let (C, W*) 
be the unique solution of (4.12)-(4.13) so that 


Then we have 
and 


(C,W*)(-5i) = (u,W)(-Ji). 

|mi(t) - Ci(t)| < Clc/llloglc/lP, T e [-( 5 i,( 52 ], 


|u(^ 2 ) - C(< 52 )| < CIc/r, |W( 52 ) - w*( 0 +)| < C|c/| |log |c 7 |r . 

Or, alternatively, assume |m(( 52 )U < 1 and the initial condition on (y,yV*) at t = 62 - Then the 
same estimate holds only at t = —Si rather that 62 - 

Proof. Denote by (V. W.,.) the solution to (4.12)-(4.13) defined in Lemma 4.4 and put C := B~^y_. 
On the one hand, from (4.13) and Corollary 4.2, we have 

w*(0-) = i(wjo+) + w,(0-)) - i(wjo+) - WJO-)) 


= W(0) - Sgn(c/) 


7ra^{72(0) 


2t/i(0) ’ 

where := Ci(0) = D]^(0). On the other hand, (4.21) and Corollary 4.2 imply 

,7rf72Ul,„,, ^ I ||3 


|W(-(5i) - >V(0) + sgn (c/)- 


2Ui 


< C\ci\ |log|c 7 | 


Furthermore, using the estimate on ui — Cj given in Lemma 4.4, we see that 
IW 40 -) ->V(0) +sgn(c 7 )^!|^( 0 )| < C\ci\ |log|c 7 ||^. 
Therefore, combining these observations, one arrives at 

|W,(-(5i) - W,(-5i)| < C|C7| |log|c7|^ 

and thus the corollary follows from Lemma 4.4 and Corollary 4.2. 


□ 


Changing from the r variable back to the X 2 variable, we obtain the other form of the limit 
system for (u*, W*)(x 2 ) = (D, W*)(T(a: 2 , s)) on X 2 € [so - <5, sq + (5]\{so}: 


(4.22) 


uC = 2 u*2 


^*2 


= k' 


w; 


17+ — Re ( 


■^^*3 — 2 ( + 


U'l 


[/+ — Re c 


+ U*3 

U*2 


w: = 0 
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where we recall Rec = U+{s). To show that the correct conditions hold at X 2 = s, we note that for 

141 < 1 

t{s±X 2 ,s) = 17+(s±4) - U+{s) = ±C/+(s)x 2 +0(141^): 
so u^, := (u*i,u* 2 ,R* 3 )^ satisfies 

M*(s + X 2 ) — u^{s — x' 2 ) = V (t(s + X 2 , s)) — V(t(s — X 2 ,s)) 

= V{t{s + X 2 ,s)) -V{t{s-X 2 ,s)) + 0(141 log^ 141)- 


Then, since, 

\Vr\ = 0(|r|-i |log|r| 
we may conclude 


T(s-4,s) + 'r(s + 4>s) =0(141 ) =0(4), 


lim (u*i,2(s + 4 ) ~'“*1,2(3 — X2)) = 0 , 


(4.23) 


W*(s+) - W*(s-) = sgn(c/) 


7 r[/"(s)u*i(s) 

|o;(3)i ^ 


irU^lis) 


lim (m*3(s + X 2 ) - m*3(s - 4)) = sgii(cj) ,| (W*(3+) + W*(s-)). 

I ®2^o+' |0+(s)r 


At last, we can state the following result on the convergence of solutions of (4.6) to those of the 
limiting system (4.22)-(4.23). This is simply a rephrasing of the previous corollary in light of the 
change of variable computations above. 


Proposition 4.2. Given any solution (u, W) of (4.6) on [s — i5, s + d] such that |u(s + d)]™ < 1, 
there exists a unique solution (u*, W*) of (4-22-4-23) so that 

(u*, W*)(s + 6) = {u, >V)(s + (5). 

In addition we have 


|ui(t) - w*i(r)| < C\ci\\ log|c/|4 T G [s - (5, s + d], 

and 

|u(s -6)- M*(s - (5)1 < C\ci\°‘, |>V(s - S) - >V*(s-)| < C\ci\ |log|c/||^ . 

Alternatively, we may assume |u(s + (5)|u; < 1 and impose the initial condition on (u*,VV*) at 
X 2 = s — 5. Then, one has the estimates above at X 2 = s -\- 5 rather than X 2 = s — S. 


4.4. Linear instabilities due to critical layers. Suppose c* G R is a regular value of C/+ on 
[0,/i+] satisfying (4.1). As in the statement of Theorem 4.1, let 

17+^({c*}) =: {si < ... < Sm} C (0, h+), where t/+(sj) ^ 0. 

We consider solutions to (4.6) for c = cr + icj G C\M close to c*. Eventually we restrict our 
attention to c* = Cfc, but the following analysis holds for any regular value c*. 

Proposition 4.3. For any a G (0, 1) and fc G N, there exist C, eo, 5 G (0, 1] depending only on a, 
k, |17+|c3([o,/(+]); and maxj |?7)|_(sj)|“^, such that 

0<si— 6 < Si + 6 <S 2 — S<...<Sm — d<Sm + d<h+ 

and the following estimates hold for any c = CR + icj satisfying |c —c*|,|c/| G (0,eo)- LG (w,W) be 
the solution of (4.6) with the initial condition 

n(h+) = (0,0,1)^, W(h+)=0 
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and let (ii», W*) 

be the solution of (4.22) for X 2 ^ 17+ 

■'({cr}) 

satisfying 


u 

*(/!+) = (0,0,1)^, 

m{h+) 

= 0 

and (4.23) at all 

s G C/+^({cii}). 

Then we have 



(4.24) 

Wi{- 

r) - m*(T)| < C\ci[ 

\ tG 

[0,h+], 

and 





(4.25) 

11 — It* + w — 

m\<C\ci\°‘, on 

[0,h+]\ 

TYl 

\J{sj-6,Sj+6) 


It is clear that the above solution {u,W) corresponds to a solution of (3.4). Moreover the 
existence and uniqueness of (u*, W*) is due to Corollary 4.2. 

Proof. When 0 < e <C 1, for any c = cn + icj in an e-neighbor hood of c*, there exist < ... < 

G (0, h+) close to Si < ... < Sm (with distance of the order of 0{e)) such that 

(4.26) {s'l,..., s:„} = {s G [0, h+] I U+{s) = cr}, C/;(s') = [/+(s,) + 0{e) ^ 0. 

We first fix (5 > 0 and take Cq sufficiently small so that |s' — Sj| < 5/2 and the estimates in 

Proposition 4.2 hold for s = s', 1 < j < to. 

For notational convenience, let sq := —5 and Sm+i '■= + 5. By induction on j from j = m 

to j = 0, we will prove the desired estimates (4.25) for u — and W — W* on the interval 
[sj + 6, Sj+i — 5] and the estimate (4.24) for m — ui* on [sj + 5, s^+i + 5] fl [0, h+]. Together, these 
imply the proposition. 

Consider j = to. Since (u, >V)(/i+) = (u*, >V*)(/i+) and (4.7) is a regular perturbation of (4.22) 
on [sm + 5, h+], as e —>• 0+, inequalities (4.25) and (4.24) hold on [sm + 5, h+] automatically, and 
actually do so with a better bound C|c/|. 

Suppose we have proved (4.25) on [sj' + 5, s^'+i — 5] and (4.24) on [sj' + S,Sj'+i + 5] fl [0,/i+], 
for j' = J + 1,..., TO and 0 < j < to — 1. We will now prove them for j' = j. Let (u, W) be the 
solution of (4.22) with (ft,yV)(sj-|-i + 5) = (it,yV)(sj-|-i + 5). On the one hand, from Proposition 
4.2 and the induction hypothesis, we first obtain 

1(11 - ii)(s,+i - 5)1 + |(W - W)(s,+i - 5)1 < C|czr, 

Km -Mi)(r)| < C\ci\“, T G [sj +1 -5, Sj+i +5]. 

On the other hand. Corollary 4.2 and the induction hypothesis imply 

|(it* - i()(sj+i - 5)1 + |(W* - 'bV)('Sj+i - 5)1 < Clc/I". 

Moreover, the boundedness of B{t) established in Lemma 4.2, along with Corollary 4.2, yields 
|(mi - Mi)(r)| < C\ci\°‘, T G [Sj+i - 5,Sj+i + 5]. 

Therefore (4.25) holds at X 2 = Sj+i — 5 and (4.24) holds on [sj+i — 5, Sj+i + 5]. Finally, since (4.6) 
is a regular perturbation of (4.22) on [sj +5, Sj+i — 5] as £ —;► 0+, inequalities (4.25) and (4.24) hold 
on [sj + 5, Sj+i — 5] and thus we obtain (4.25) on [sj + 5, Sj+i — 5] and (4.24) on [s^ + 5, + 5]. 

The proof of the proposition is then completed by induction. □ 

In order to obtain a solution of (3.4) and (3.5), we need to prove iti(O) > 0 which follows from 
Proposition 4.3 and the next lemma. 
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Lemma 4.5. Assume C/+ € Let J7_,_^({cfl}) = and say that is the 

solution of (4.22) on [0, h+\ \ {s(,..., s'^} that satisfies 

u4h+) = {0,0,lf, W4h+) = 0 

and (4.23) at any s € [/^^({cij}). Then (m*(0), W*(0)) is in cr for cr in a neighborhood of c^. 
Moreover 


M*i,3(2^2) > 0, a:2 G [0,/i+]; u*i(s(„) > 0, and if m> 2, then u^:i{s'^_i) > 0 

where s'j, j = 1,... ,m, are defined in (4.26). 


Proof. Since the transformation matrix B{t) given in Lemma 4.2 is smooth in both r and the pa¬ 
rameter s close to So (in terms of the notations in Lemma 4.2), Corollary 4.2 implies the smoothness 
in cr of (u*, W*) near (but not at) a singularity. Equation (4.22) is regular away from {s(,..., sl^} 
with uniform bounds on the coefficients. Therefore we see that (u*(0), W(0)) is in cr near c*. 

To prove > 0, notice Lemma 4.3 and the boundary conditions for (u*, W*) at h+ imply 


(4.27) 


~ =0 on [ 0 , h+]. 


It follows immediately that and u *3 can not vanish simultaneously at any X 2 ^ C/l|r^({c/j}), and 
this occurs only when the solution is trivial. Conservation law (4.27) actually further implies that, 
if one of 21 * 1 ^ 3 ( 0 : 2 ) = 0 at some X 2 ^ C7ip^({cfl}), it does not change sign since the other one does 
not vanish in a neighborhood. Moreover, since is continuous on [0,h+], if u*i(s') ^ 0, then 
u *3 does not change sign near this s' due to (4.27). Finally, if m*i(s') = 0, Lemma 4.3 implies 
(u*,yV*) is smooth near this s' and u* 3 (s') ^ 0. Therefore it*i does not change sign near s' as 
well. Summarizing the above discussion, we conclude that u*i ,3 > 0 on [0,/i+J. 

Finally we prove m*i(s(„) > 0 and u*i(s(„_i) > 0 if m > 2. Even though we will continue to work 
in the framework of (4.22), the calculation is essentially carried out to the form of the equation 
used in the proof of Lemma 3.1 where tp = ?//([/+ — c) was considered. Let 

U' 

H{x2 ) ■■= M ,2 - 77 - X2 G [0, /l+]\[/r^({Cfl}). 

C+ — CR 


On the one hand, one may compute 

H'{x2) = -b ( 


U+ - Cr 

> -b ( 


c; .2 2c/; 

j -b u*3 - 77 -U *2 


u+ - Cr 


U+ - Cr 

ifl- > 0 , 


in view of (4.27). Indeed, the conservation law (4.27) also implies that the above derivative vanishes 
only for the trivial solution. Since H{h+) = 0, the monotonicity of H implies 


(4.28) 


lim H{x 2 ) G [— 00 , 0 ). 

X2^s'^ + 


On the other hand, from (4.27), it is clear H{x 2 ) = 0 if m*i(x 2 ) = 0 at some X 2 G [0, /i+]\[/i|r^({c_R}). 
Moreover, suppose m*i(s') = 0 for some j = 1,..., m. Lemma 4.3 and (4.22) imply M*(a: 2 ), and thus 
H{x 2 ) as well, is smooth near X 2 = s'y Moreover (4.27) yields 11 * 2 ( 5 )) = W*(s'±) = 0. Therefore 
^* 1 ( 5 )) = 2u*2(s') = 0 which leads to H{s'j) = 0. Consequently, (4.28) implies m*i(s(*j) > 0. From 
(4.23), we obtain >V*(s(„—) ^ W*(sm+) = 0. Again, Lemma 4.3 implies u*i(s(*j_i) 7 b 0 if m > 2. 
This completes the proof of the lemma. □ 
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Corollary 4.4. Assume in addition to the hypotheses of the previous lemma that Uf_{sj), j = 
1,... ,m, are all non-positive or all non-negative, and U'f_{sjg) 0 for jo = m or m — 1. Then, for 
any cr near c*, there exists a unique solution of of (4.22) such that 

u*i(0) = 1, m*i,2(/i+) = 0, W,(/i+) = 0, 

and (u*, yV*)(0) is in its dependence on cr. Furthermore, m*i(x2) ^0 at X 2 = 0, s',..., and 


W*(0) 


sgn (C7)7r^ 
j=i 


U!^(s'Ms') 


^0. 


The additional assumption and Lemma 4.5 imply W* 0 and W* does not change sign on 
2^2 < Consequently, (4.27) and Lemma 4.5 imply u*i(x 2 ) > 0 for < s'^^. The corollary 
follows from normalizing the solution furnished by Lemma 4.5. 

Finally, combining Proposition 4.3 and Lemma 4.5 gives the following statement about the 
solution of Rayleigh’s equation (3.4) with a near singular coefficient. 


Proposition 4.4. Under the assumptions in Corollary for c = cr -\- icj sufficiently close to 
c^, with c/ 7 ^ 0, there exists a unique solution y to (3.4) and (3.5), which corresponds to a unique 
solution (u, W) of (4.6) satisfying 


(4.29) wi(0) = l, ui,2(/i+) = W(/i+) =0, U3{h+)>0. 

Let (it*,yV*) be the solution given in Corollary for cr and sgn(c/). For any a € (0,1), there 
exists C > 0 such that 

\u-u4 + |W-W,| < C|c/|“. 

Moreover 

y'(0) = u2(0)+iW(0). 

Proof. According to (4.4), any solution y of (3.4) clearly gives rise to a solution {u,yV) of (4.6) 
with the above specihed boundary conditions. Such a solution (u, W) is unique, much as we saw in 
Corollary 4.4, in light of Proposition 4.3. 

Conversely, we can reconstruct y from such (u,W). In fact, the conservation law (4.5) implies 
Ml and M 3 can not vanish simultaneously, unless the solution is trivial. This and (4.5) further imply 
that Ml ,3 do not change sign on [0, h+] and thus both remain nonnegative. Let 9 = 9{x2) be defined 
as the solution of 

O' = —, 6>(0) = 0. 

Ml 

This function is well-defined as, if Mi(a; 2 ) = 0 at some X 2 £ [0, h+], then (4.6) and conservation law 
(4.5) imply 

Ui(a^ 2 ) = 2 m 2 (x 2 ) = 0 = W(x 2 ) = W'(a; 2 ), and m'/(x 2 ) = M 3 (x 2 ) > 0. 

Let y := y'Mie'®. One may compute using the definition of 9, (4.5), and (4.6), that 

y' = ^L(m2 +iW)e^^, 
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while 


y” = (- M2M1 " (m 2 + *>v) + + U3 + 


icjU" 


{U+ - Cr) + cj 


{u+ - Cr) + C 


:U1, 


+ iWu^ ^ {u2 + iW )) e*' 
1 


= ( - ujHui + W^) + + »3 + J7^)“.))«" = ('•■" + 

Therefore y solves (3.4). The estimates on y are from Proposition 4.3. 


□ 


Proof of Proposition 4.1. The correspondence between the solution y of (3.4) and (u, W) was 
established in the proof of Proposition 4.4,. The properties of (u, W) and the convergence estimates 
were already obtained in the previous lemmas and propositions. To complete the proof, one only 
needs to confirm that the jump conditions on ( 2 /*,y() are satisfied at each s where U+{s) = c*. In 
fact, the log^ growth bound on |m» 3 | implies y'^ has at most a logarithmic singularity, which shows 
that y* is Holder continuous. This allows us to infer that in the limit X 2 —t 0+, 

{y'*y){s + X 2 ) - {yiy){s - X 2 ) = {y'^{s + X 2 ) - yi{s - X 2 ))y*{s - X 2 ) 

- y'(s + X 2 ){y*{s + X 2 ) -y*{s- X 2 )) 

-^y*{s) lim (y'(s + a; 2 ) - y'(s - X 2 )). 

X 2 —k 0 + 

As and W* are the real and imaginary parts of y(,y*, condition (4.2) follows from the above 
calculation and (4.23). □ 


Finally we are in the position to prove the main theorem. 


Proof of Theorem 4.1. Let {u^,yV^) be the solution of the limiting system (4.22) with the 
parameter cr = Ck and put 


(4.30) 


(t/+(0)-Cfc)^ - t/"(s,)u#i(s,-) 

'^2cfe|fc|tanh(|fc|/i_) ^ |C^+(®i)l 


Note that > 0 due to Lemma 4.5. Define a mapping G = (Gi, G 2 )(ci, C 2 , e) by 


Gi(ci,C 2 ,e) := Cfc + ci - /i^(eRey'(O), elmy'(0),e) 
G 2 (ci,C 2 ,e) := c# + C 2 - Imy'(0)//(eRey'(0), elmy'(O), e) 


where /rj are given in (3.8) and y is the solution of (3.4) and (3.5) with the parameter 
(4.31) c = Cfc + Cl + ie(c^ + C 2 ). 

The existence and uniqueness of y is ensured by Proposition 4.4 for small ci and e. Clearly G is 
smooth for e(c^ + £ 2 ) > 0. In addition, the zero-set of G corresponds to the solutions of (3.8), and 
thus solutions to (3.4)-(3.6). Proposition 4.3 and Corollary 4.4 imply 

y'(0) = U 2 ( 0 ) -h tW(0) = m* 2 ( 0 ) ^W*(0) 0(e“), 

where a can be taken arbitrarily in (0,1), {s)^,..., := U^^{{ck + ci}), and (u*,W*) is the 

solution of (4.22) and (4.23) with the parameter Cfc -I- Ci. From the smoothness of (m*(0), W*(0)) 
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in Cl due to Corollary 4.4, we have 

y'( 0 ) =u^ 2 ( 0 ) + iW#( 0 ) + 0 (|ci| + e“) = m# 2 ( 0 ) - itt^ 

On the other hand, (3.8) and (3.9) give us that 

G(ci, C 2 , e) = (ci + 0 (e), £2 + 0 (e“ + |ci|)). 

Take (£ 1 , 62 ) in a rectangle [—15, d] x [—M(5, Md], where the fixed constants 5 is small and M is 
large. By considering the the image of the boundary of this rectangle under the mapping G(-, •, e), 
a standard degree theory argument implies that G has a zero point near 0 for any fixed small e > 0. 
Unravelling definitions, this corresponds to an eigenvalue near Cfc with positive imaginary part. The 
proof of the theorem is complete. □ 


Remark 4.2. 1.) To ensure the existence of an instability, the above argument makes clear that one 
needs only that 


(4.32) 



ic^;(si)i 


> 0 , 


and thus those given in Theorem 4.1 are sufficient, but not necessary.. 

2.) If both the positive and negative values of belong to the range of C/+ and satisfy the 
assumption in Theorem 4.1, the above proof implies that there exist at least two distinct unstable 
modes. 


5. Instability by other means 

We can summarize the central conclusion of the previous two sections as follows. Assume t/+ S 
and fix a wave speed k. Let Ck denote the corresponding wave speed for the capillary-gravity 
water wave beneath vacuum problem (i.e. e = 0) given in (3.7). On the one hand, if a sequence of 
unstable wave speed Ck,e„ of the water-air problem converges to Ck as the density ratio e„ O-f, 
then Lemma 3.1 implies that there must be a critical layer in the shear flow in the air, i.e., Ck € 
i7+([0, h+]). On the other hand, under the non-degenerate shear condition on U\. and some sign 
condition on [/", unstable wave speeds may bifurcate from Ck at e = 0. 

The situation for the Kelvin-Helmholtz instability is different. Inequality (3.11) indicates that 
Kelvin-Helmholtz instability occurs if and only if the parameters (fc, e) are in a region K C {fc, e > 
0}. With surface tension, there is a positive distance between K and the set {e = 0}, that is, for any 
wind speed Uo > 0, instability does not occur if e is too small. Without surface tension, on the one 
hand for any fixed fc > 0, we have (fc, e) ^ K, i.e., the wave number k is stable if e > 0 is sufficiently 
small. In this case, no unstable wave speed bifurcates from Cfe at e = 0. On the other hand, the 
distance between K and {e = 0} is zero, because any sequence {(fc„, Cn)} C K with —?► 00 must 
satisfy e„ —>■ 0. One implication of this is that, for any e > 0, all sufficiently large modes k are 
unstable. In contrast to the water-vacuum setting, this instability “bifurcates from infinity.” 

Relaxing the assumption that U+ G may lead to additional instability that does not fall into 
the critical layer theory. To illustrate this point, in this section we show that there exist background 
flows U that are linearly unstable at a wave number k for which Ck is in the range of 17+, but the 
critical layer and support of f/" are separated by a distance uniform in small e > 0. In fact, these 
are solutions of the Euler system for which U” is a 5-measure of negative mass, and the critical 
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layer is at an inflection point. In the view of Remark 2.3, it is justified to use (3.4)”(3.6) to study 
the instability of such non-smooth shear flows. 

For simplicity, in these computations we take cr = 0 and h± = oo. Consider background profiles 
of the form 


(5.1) 


U+{x2) := 


IJ,X2 0 < X2 < X2 
11 X 2 X2 > X2- 


Here ^ and X 2 are parameters that we may choose freely. Notice in particular that 


Clearly the value 

will play an important role. 


u'l = 




[/* := U+{x*2) = fiX2 


Proposition 5.1. Fix a wave speed k. Let a = 0, and take h± = 00 . There exists U € C°’^(K) 
of the form (5.1) that are unstable in the sense that, for aZ/ 0 < e <C 1, there exists {y,c) solving 
(3.4)-(3.6) with 

0 < Imc = 0{y/e). 

Moreover, for this profile, C/* — \fgjk has a positive lower bound uniform in e, hence the critical 
layer occurs at an inflection point of U away from X 2 ■ 


It is worth pointing out that the exponential growth rate 0(i/e) of the instability here is much 
greater than 0 (e) predicted by the critical layer theory. 


Proof. Taking U as in (5.1), it is possible to explicitly solve the Rayleigh equation (3.4) for y. In 
general, we find 


(5.2) 


r Aie-'^^^ X2 > x*2 

\ A2e-^^^ + Ase'^^-^ 0 < X2 < x*2 


where Ai, A 2 , H 3 are given by the following linear system: 


A 2 -h H 3 — 1 

A 1 -A 2 - Ase^^^^ = 0 

- k)Ai + kA2 - kAse^^^'^ = 0. 

— c 


The first of these is from the boundary condition at X 2 = 0, the second is to ensure continuity 
across and the third arises from the jump condition for y' across X 2 . Solving the system yields 

^ __ 2fc(17, - c) _ 

^ 2 fc([/* — c) + /i(e “^^'^2 — 1 ) 

2 k{U, - c) - M 

^ 2fc(t/* — c) + — 1) ’ ^ 2fe(C/* — c)-I-— 1) 


Inserting these computations into (5.2), we see that 


y'{0) = k{A3-A2) = -k^, 

c-P 

where 

(5.3) a:=C/*-^(l + e-2^-^), /? = I/, - A(i _ ). 
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With this in hand, the dispersion relation (3.6) becomes 


g{l -e) = ekc^ 


Note that, by (5.3), 


= C/, 


1 - 


c — 

1 

2kxi 


+ (?k — ec/r. 

(l_e-2fc^5) . 


The quantity in square brackets on the right-hand side above has range (0,1), with 1 and 0 being 
the limits as is taken to -l-oo and 0-1-, respectively. It is therefore easy to see that there exist 
many choices of the parameters X 2 and /i for which 


/3 = 


< U, 


and 


/r > 0. 


The dispersion relation can then be simplihed into the following polynomial: 


(5.4) 


0 = /(c,e) := c- 


\/f)( 




9 


c"-efc-f(l-e) +ec"(c-a). 


Notice that c = \fgjk is a double root of /(•, 0). Indeed, 


/(l/p0) = (^c/)( 


0 )= 0 , ( 9 . V )(. ^, 0 )= 4 , ^>0 


(5./)( 


!■“>=!< 


f -»)=> 0^ 


Since / is a polynomial of c, these facts imply that, for 0 < e <C 1, there exist a complex conjugate 
pair c± of solutions to f{c, e) = 0 with 


Rec± — + 0(e), 


0 < ±Imc± = O(-ye). 


□ 
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